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Abstract 

Recently we have developed a novel chiral power counting scheme for an eflective field theory of nuclear matter 
with nucleons and pions as degrees of freedom [1] . It allows for a systematic expansion taking into account both 
local as well as pion-mediated multi-nucleon interactions. We apply this power counting in the present study to 
the evaluation of the pion self-energy and the energy density in nuclear and neutron matter at next-to-leading 
order. To implement this power counting in actual calculations we develop here a non-perturbative method 
based on Unitary Chiral Perturbation Theory for performing the required resummations. We show explicitly 
that the contributions to the pion self-energy with in-medium nucleon-nucleon interactions to this order cancel. 
The main trends for the energy density of symmetric nuclear and neutron matter are already reproduced at 
next-to-leading order. In addition, an accurate description of the neutron matter equation of state, as compared 
with sophisticated many-body calculations, is obtained by varying only slightly a subtraction constant around 
its expected value. The case of symmetric nuclear matter requires the introduction of an additional fine-tuned 
subtraction constant, parameterizing the effects from higher order contributions. With that, the empirical 
saturation point and the nuclear matter incompressiblity are well reproduced while the energy per nucleon as a 
function of density closely agrees with sophisticated calculations in the literature. 



1 Introduction 



In the last decades Effective Field Theory (EFT) has been applied to an increasingly wider range of phenomena, 
e.g. in condensed matter, nuclear and particle physics. An EFT is based on a power counting that establishes a 
hierarchy between the infinite amount of contributions. At a given order in the expansion only a finite amount 
of them has to be considered. The others are suppressed and constitute higher order contributions. In this work 
we employ Chiral Perturbation Theory (CHPT) [2 4] , which is the low-energy EFT of QCD and takes pions (and 
nuclcons as well for our present interests) as the degrees of freedom. CHPT is related to the underlying theory of 
strong interactions, QCD, because it shares the same symmetries, their breaking and low-energy spectrum. It has 
been successfully applied to the lightest nuclear systems of two, three and four nucleons [5-11]. Nonetheless, still 
some issues are raised concerning the full consistency of the approach and variations of the power counting have 
been suggested [12-19]. A common technique for heavier nuclei is to employ the chiral nucleon-nucleon potential 
delivered by CHPT in standard many-body algorithms [20,21], sometimes supplied with renormalization group 
techniques [22, 23]. One issue of foremost present interest is the role of multi-nucleon interactions involving three 
or more nucleons in nuclear matter and nuclei [11,20,23-25]. 

Ref. [26] derived many-body field theory from quantum field theory by considering nuclear matter as a continuous 
set of free nuclcons at asymptotic times. The generating functional of CHPT in the presence of external sources 
was deduced, similarly as in the pion and pion-nucleon sectors [27,28]. These results were applied in ref. [29] to 
study CHPT in nuclear matter, but including only nucleon interactions due to pion exchanges. Thus, the local 
nucleon-nucleon (and multi-nucleon) interactions were neglected. This approach was later extended in ref. [30] 
to finite nuclei and the pion- nucleus optical potential is calculated up to 0{p^). In ref. [1] an extended power 
counting is derived that takes into account simultaneously short- and long-range multi-nucleon interactions. Notice 
that many present applications of CHPT to nuclei and nuclear matter [25,29,31-39] only consider mcson-baryon 
Lagrangians. Short-range interactions are included without being fixed from the free nucleon-nucleon scattering. 
E.g. [35] fits the in-medium local nucleon-nucleon interaction, in terms of just one free parameter, to reproduce 
the saturation properties of symmetric nuclear matter. In addition, the nucleon propagators do not always count 
as 1/fc, with k a typical nucleon three-momentum, but often they do as the inverse of a nucleon kinetic energy, 
m/fc^ (with m the nucleon mass), so that they are unnaturally large. This is well known since the seminal papers 
of Weinberg [3,4]. This fact invalidates the straightforward application of the pion-nucleon power counting valid 
in vacuum as applied e.g. in refs. [29,31-33,40]. 

We implement here non-perturbative methods to perform actual calculations employing the power counting of 
ref. [1], which requires the resummation of some scries of in- medium two- nucleon reducible diagrams. We employ 
the techniques of Unitary CHPT (UCHPT) [41-44] that are extended to the nuclear medium systematically in 
a way consistent with the chiral power counting of ref. [1]. Our theory is applied to the problem of calculating 
up to ncxt-to-lcading order (NLO) the pion self-energy and energy density in asymmetric miclcar matter. The 
former problem is related to that of pionic atoms since the pion self-energy and the pion-nucleus optical potential 
are tightly connected [45,46]. The issues of the pion-nucleus S-wave missing repulsion, the renormalization of the 
isovector scattering length a~ in the medium [37, 47] and the energy dependence of the isovector amplitude [46] are 
not settled yet, despite the recent progresses [29,46,48]. Ref. [1] found that the leading corrections to the linear 
density approach for calculating the pion self-energy in nuclear matter are zero. We show here these cancellations 
explicitly within the developed non-perturbative techniques. We also show the related cancellation between some 
next-to-next-to-leading order (N^LO) pieces. The calculation of the energy density of nuclear matter starting 
from nuclear forces is a venerable problem in nuclear physics [31,33,35,49-54]. Our calculation of the energy 
density £ to NLO already leads to saturation for symmetric nuclear matter and repulsion for neutron matter. 
Indeed, an accurate reproduction of the equation of state for neutron matter can be achieved by varying slightly 
one subtraction constant around its expected value. For the case of symmetric nuclear matter an additional 
fine-tuning of a subtraction constant is necessary to obtain a remarkable good agreement between our results 
and previous existing sophisticated many-body calculations [49,54,55]. The saturation point and nuclear matter 
incompressibility are reproduced in good agreement with experiment. 

After this introduction, we briefly review in section 2 the novel chiral power counting in the medium developed 
in ref. [1]. The contributions to the pion self-energy in the nuclear medium that arise from tree- level pion-nucleon 
scattering diagrams and from the one-pion loop nucleon self-energy are the subject of section 3. We dedicate 
section 5 to the evaluation of the part of the pion self-energy due to the dressing of the nucleon propagators in 
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the medium because of the nucleon-nucleon interactions. For their calculation one requires the nucleon-nucleon 
scattering amplitudes in the nuclear medium which are calculated in the preceding section 4. The terms of the 
pion sclf-cnorgy chic to the nucleon-nucleon interactions that are not part of the nucleon self-energy are calculated 
in section 6 at NLO, where we also give some N^LO contributions. The derivation of the necessary loops involved 
in this calculation is performed in Appendix B. Section 7 is dedicated to the evaluation up to NLO of the energy 
density. Section 8 contains a short summary and the conclusions. In the Appendices we derive various results that 
are used in the main body of the paper. Appendix A offers a derivation of the partial wave expansion of nucleon- 
nucleon scattering in the nuclear medium and vacuum. The Appendices C, D and E develop the calculation of the 
in-medium integrals needed for the evaluations performed in other sections. 



2 In-medium chiral power counting 

We briefly review the chiral power counting for nuclear matter developed in ref. [1]. Let us start by introducing 
the concept of an "in-medium generalized vertex" (IGV) given in ref. [26]. Such type of vertices arises because one 
can connect several bilinear vacuum vertices through the exchange of baryon propagators with the flow through 
the loop of one unit of baryon number, contribiited by the nucleon Fermi-scas. At least one Fermi-sca insertion is 
needed because otherwise we would have a closed vacuum nucleon loop that in a low-energy effective field theory is 
completely decoupled. It is also stressed in ref. [29] that within a nuclear environment a nucleon propagator could 
have a "standard" or "non-standard" chiral counting. To sec this note that a soft momentum Q ^ p, related to pions 
or external sources can be associated to any of the vertices. Denoting by k the on-shell four-momentum associated 
with one Fermi-sea insertion in the IGV, the four-momentum running through the j*'* nucleon propagator can be 
written as pj = k + Qj, so that 



{k + Qjf -w?+ie Q] + 2Q0.E{k) - 2Q,k + ie ' 

where i?(k) = k^/2m, with m the physical nucleon mass (not the bare one), and is the temporal component of 
Qj . We have just shown in the previous equation the free part of an in-medium nucleon propagator because this is 
enough for our present discussion. Two different situations occur depending on the value of Q^. If Qj = 0{m^) = 
0{p) one has the standard counting so that the chiral expansion of the propagator in eq. (2.1) is 

Thus, the baryon propagator counts as a quantity of 0{p^^). But it could also occur that Q° is of the order of a 
kinetic nucleon energy in the nuclear medium or that it even vanishes.*^ The dominant term in eq. (2.1) is then 



' Q2 2Qj • k - ie 



-^ o2 1, .V ' (2-3) 



and the nucleon propagator should be counted as 0(^3"^), instead of the previous 0(p~^). This is referred to as 
the "non-standard" case in ref. [29]. We should stress that this situation also occurs already in the vacuum when 
considering the two-nucleon reducible diagrams in nucleon-nucleon scattering. This is indeed the reason advocated 
in ref. [3] for solving a Lippmann-Schwinger equation with the nucleon-nucleon potential given by the two-nucleon 
irreducible diagrams. 

In order to treat chiral Lagrangians with an arbitrary number of baryon fields (bilinear, quartic, etc) ref. [1] 
considered firstly bilinear vertices like in refs. [26, 29], but now the additional exchanges of heavy meson fields of 
any type are allowed. The latter should be considered as merely auxiliary fields that allow one to find a tractable 
representation of the multi-nucleon interactions that result when the masses of the heavy mesons tend to infinity. 
Note that such methods are also used in the so-called nuclear lattice simulations, see e.g. [56]. These heavy meson 
fields are denoted in the following by and a heavy meson propagator is counted as 0(p°) due to the large meson 



An explicit example is shown in section 6 of ref. [29] 



3 



mass. On the other hand, ref. [1] takes the non-standard counting case from the start and any nucleon propagator 
is considered as 0{p~^). In this way, no diagram, whose chiral order is actually lower than expected if the nucleon 
propagators were counted assuming the standard rules, is lost. In the following ^ kp ^ C'(p) ^-re taken of 
the same chiral order, and are considered much smaller than a hadronic scale A^^ of several hundreds of MeV that 
results by integrating out all other particle types, including nucleons with larger three-momentum, heavy mesons 
and nucleon and delta resonances [4] . The formula obtained in ref. [1] for the chiral order of a given diagram is 

V m Vp 

v = A-E + Y,{n^+^^-^)+Y.^d,+u:,-l)+Y.^v,-l) + Y.''i ■ (2.4) 

i=l i=l i=l i=l 

where E is the number of external pion lines, ti; is the number of pion lines attached to a vertex without baryons, 
^i is the chiral order of the latter with 14 its total number. In addition, di is the chiral order of the i*'' vertex 
bilinear in the baryonic fields, Vi is the number of mesonic lines attached to it, that of only the heavy lines, V is 
the total number of bilinear vertices, Vp is the number of IGVs and m is the total number of baryon propagators 
minus Vp, V = Vp + m. The previous equation can be also written as 

v„ V 

u = A-E + Yl{ni + - 4) + Y^{di +Vi+uJi-2) + Vp . (2.5) 

i=l i=l 

It is important to stress that v is bounded from below as explicitly shown in ref. [1] . Because of the last term in 
eq. (2.4) adding a new IGV to a connected diagram increases the counting at least by one unit because Vi > 1. The 
number v given in eq. (2.4) represents a lower bound for the actual chiral power of a diagram, jjb, so that n> u. 
The actual chiral order of a diagram might be higher than v because the nucleon propagators are counted always 
as 0{p~'^) in eq. (2.4), while for some diagrams there could be propagators that follow the standard counting. 
Eq. (2.4) implies the following conditions for augmenting the number of lines in a diagram without increasing the 
chiral power by adding i) pionic lines attached to lowest order mesonic vertices, £i = rii = 2, ii) pionic lines attached 
to lowest order meson-baryon vertices, di = Vi = 1 and iii) heavy mesonic lines attached to lowest order bilinear 
vertices, di = 0, = 1. One major difference between our counting, eq. (2.5), and Weinberg one [3,4] is that ours 
applies directly to the physical amplitudes while the latter applies only to the potential. 

We apply eq. (2.4) by increasing step by step Vp up to the order considered. For each Vp then we look for those 
diagrams that do not further increase the order according to the rules i)-iii). Some of these diagrams are indeed 
of higher order and one can refrain from calculating them by establishing which of the nucleon propagators scale 
as 0{p~^). In this way, the actual chiral order of the diagrams is determined and one can select those diagrams 
that correspond to the precision required. 

It is worth realizing that eq. (2.5) can be also applied in vacuum in order to determine the relative weight of the 
different diagrams. In this case, the needed Fermi-sea insertion for each IGV is split in two external nucleon lines, 
both in- and out-going ones. For vacuum Vp is constant because in the EFT there is no explicit closed nucleon 
loops and baryon number is conserved. As stressed above, the expressions between brackets in eq. (2.5) do not 
increase despite the diagrams become increasingly complicated. As a result, one takes Vp constant and determines 
the leading, next-to-leading, etc, contributions as indicated in the previous paragraph. In order to derive eq. (2.5) 
in ref. [1], a term 3Vp was summed because for each IGV there is a Fermi-sea insertion. Since in vacuum there is 
no sum over nucleons in the sea, one should subtract this contribution so that we would have 4 — 2Vp instead of 
4: + Vp. However, this just modifies the absolute order of a diagram but not the relative one between contributions 
which remains invariant, and this is what matters for explicit calculations. 



3 Meson-baryon contributions to the pion self-energy 

Here we start the application of the chiral counting in eq. (2.4) to calculate the pion self-energy in the nuclear 
medium up to NLO or 0{p^). The different contributions are denoted by H, and are depicted in fig. 1. In terms 
of the pion self-energy H the dressed pion propagators reads 
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Figure 1: Contributions to the in-medium pion self-energy Hi up to NLO or 0(p"'). The pions are indicated by the dashed 
lines and the squares correspond to NLO pion-nucleon vortices. A wiggly line is the rmclcon-rmclcon interaction kernel, 
given below in fig. 5, which is iterated as indicated by the ellipsis. The thick lines correspond to closed Fermi-sea insertions, 
while the thin lines represent in-medium nucleon propagators, eq. (3.4). The external pion lines in diagrams 3a, 3b, 8 and 
10 should be understood as leaving or entering the diagrams. 



The in-medium nucleon propagator [57], Go{k)i^, is 



^fea - iki) , 0{m-u 



A;0-£;(k)-ie k'^ - E{]i.) + ie 

^ - + i27r^(Ci3 - |k|)5(fc° - E{k)) . (3.2) 



In this equation the subscript is refers to the third component of isospin of the nucleon, with = +1/2 for 
the proton and —1/2 for the neutron, and is the corresponding Fermi momentum. We consider that isospin 
symmetry is conserved so that all the nucleon and pion masses arc equal. The first term on the right hand side 
(r.h.s.) of the first line of eq. (3.2) is the so-called hole contribution and the last term is the particle part. In 
the second line, the first term is the free-space part of the in-medium nucleon propagator and the last term is 
the density-dependent one (or a Fermi-sea insertion). The proton and neutron propagators can be combined in a 
common expression 

Go{k) = (^o{^, - |k|) + ^e{^n - |k|)^ ^ 



2 -vsp . 2 ' ' V fcO - E(k) - ie 

1 

fcO - E{]s.) + ie ' 



+ ( i^e(|k| - e,) + V^^(|k| - U) ,0 J'n.^ , .V ' (3-3) 
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or in the equivalent form, 
Go(fc) 



ko - E{]s.) + ie 



1+T3 



0{^P - |k|) + 



T3 



(3.4) 



In the fohowing, cr* and r' correspond to the Pauh matrices in the spin and isospin spaces, respectively. 

For the evaluation of the diagrams 1-6 we employ the 0{p) and Heavy Baryon CHPT (HBCHPT) 

Lagrangians [58, 59] 



,(2) 
-■kN 



N 



\N ■ 



(3.5) 



where the ellipses represent terms that are not needed here. In this equation, N is the two component field of 
the micleons, gA is the axial- vector pion-miclcon coiipling and = 9^ + the covariant chiral derivative, with 
T/j — The pion fields ti{x) enter in the matrix u = exp(iT • 7r/2/), in terms of which = i 

and U = u^, with / the weak pion decay constant in the; SU(2) chiral limit. The Ci are chiral low-energy constants 
whose values are fitted from phenomenology [58] . The in-medium pion self-energy from C^n in HBCHPT has been 
calculated up to two loops in refs. [40,60]. 

The diagrams 1-6 were calculated in ref. [1]. We give here more details in their derivation. The first diagram in 
fig. 1 corresponds to 



ni = -^j2-^ij3{Pp - Pn 



(3.6) 



and arises by closing the Weinberg- Tomozawa term (WT) in pion-nucleon scattering. In the previous equation 
Pp{n) — Cp(„)/37''^ is the proton(neutron) density. Hi is an S-wave isovector self-energy. 

The diagram 2a in fig. 1 is represented by Il2a and is obtained by closing the nucleon pole terms in pion-nucleon 
scattering, with the one-pion vertex from the lowest order meson-baryon chiral Lagrangian C^j^ [58], 



2a 



9\ 



4/2 y (27r)^ 



Tr 



^(Cp-|k|) + ^^(Cn-|k|) 



tV^ u ■ q(T • q 
£;(k) - - i;(k - q) + ?:e 



(3.7) 



In eq. (3.7) we have not included the in-medium part of the intermediate nucleon propagator because g° ~ m,r ^ 
£■(!<;) — i?(k — q), so that the argument of the in-medium Dirac delta-function of eq. (3.4) cannot be fulfilled. By 
the same token 

1 _ J_ _ £;(k) - Ei^ - q) 



i;(k) - £;(k - q) - 9° 



(g0)2 



+ 0(9) , 



(3.8) 



and the 0(g) terms contribute one order higher than NLO. On the other hand, 



£;(k) - i?(k - q) 



q2 - 2k • q 

2m 



(3.9) 



and the k • q term, when included in eq. (3.7), does not contribute because of the angular integration. Then, 



Ho 



9\ 



4/2g0 



1 



d^k 



2mq° J J (27r) 



Tr 



1 + r3 



Oi^P - |k|) 



0i^n - |k|) 



r'r-' CT • qcj • q 



(3.10) 



The same procedure can be applied to the diagram 2b of fig. 1 (which corresponds to the same expression as Il2a 
but with the exchanges — >■ — and i ^ j). Summing both, one has 



^ ~ 2f2q0 ^ijSKPp - Pn) , 



2\2 



IK 



4/2 mqo 



2^ijiPp + Pn) ■ 



(3.11) 
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The superscripts iv and is refer to the isovector and isoscalar nature of the corresponding contribution to 112, 
respectively. Both are P-wave self-energies but n^'' is a recoil correction of Yi™ and it is suppressed by the inverse 

of the nucleon mass. 

The rest of the diagrams in fig. 1 are NLO. We now consider the sum of the diagrams 3a and 3b, where the 
squares indicate a NLO one-pion vertex from eq. (3.5). It should be understood that the pion lines can 

leave or enter these diagrams. We also employ the expansion of cq. (3.8) for the nucleon propagator, although for 
this case it is only necessary to keep the term ±l/g° because the diagram is already a NLO contribution. The 
calculation yields 

2 2 

^^ = ^{PP + Pn)kj . (3.12) 

This is a P-wave isoscalar contribution. In this case the NLO pion-nucleon vertex is a recoil correction of the LO 
one and this is why lis is suppressed by the inverse of the nucleon mass. The diagram 4 in fig. 1 is given by 

n4 = (2C1 ml - g2(c2 + C3 - ^) + C3 {pj> + Pn) ■ (3.13) 

114 is an isoscalar contribution in which the term —25ijCzC^ {pp+ pn) / P is P-wave and the rest is S-wave. Indeed, the 
low-energy constant C3 is known to be dominated by the contribution of the A(1232) [61]. For a Fermi momentum 
^ ~ 2m7r, corresponding to symmetric nuclear matter saturation, the Fermi energy of a two- nucleon system is 
around 80 MeV, which is still significantly smaller than the A-nucleon mass difference. One then expects that 
integrating out the A-resonance and parameterizing its effects in terms of the chiral counterterms is meaningful 
in the range of energies wc arc considering. This is indeed an important conclusion of rcf. [62] where chiral EFTs 
with/without As are employed to evaluate different orders of the two- nucleon and three-nucleon potentials. We 
leave as a future improvement of our results to include explicitly the A resonances. 

Let us consider the contributions to the pion self-energy due to the one-pion loop nucleon self-energy. This is 
represented by the diagrams 5, 6a and 6b in fig. 1. These diagrams originate by dressing the in- medium nucleon 
propagator of the diagrams 1, 2a, 2b, in order, with the one-pion loop. As a preliminary result let us first evaluate 

TT 

p ' ^ P 

k k-e k 



Figure 2: One-pion loop contribution to the nucleon self-energy in the nuclear medium. The four-momenta are indicated 
below the corresponding lines in the figure. 

the nucleon self-energy in the nuclear medium corresponding to fig. 2. First, we consider the case of a neutral 
pion. The results for the charged pion contributions follow immediately from the tt^ case. In HBCHPT the proton 
self-energy due to the one-7r° loop is given by. 



f2 - J (2Tr)D {e^-ml+ie)iv{k-£)+ie) 
+ 2n'-As,S. J ^^^-_^J(.(fc - e))Oi^^ - [k - 1[) . (3.14) 

Here, 1 is the vector made up from the spatial components of £ and v is the four-velocity normalized to unity 
(t;2 = 1)^ such that the four- momentum of a nucleon is given by p = mv + k, with k a small residual momentum 
{v ■ k <^ m). In practical calculations we take v = (1,0) and Z) ^ 4. Notice that the last integral in the previous 
equation is convergent because of the presence of the Dirac delta and Heaviside step functions. Instead of the full 
non-relativistic nucleon propagator eq (3.2), HBCHPT typically implies the so-called extreme non-relativistic limit 
in which i^(k) — )• 0, see e.g. ref. [58]. Given the properties of the covariant spin operator 5^ [58] it follows that 

= 54rr = l{{v- if - f) . (3.15) 
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For the vacuum part we then have the integral, 



4/2 J (27r)^ (^2 -rnl+ ie){v{k - £) + ie) ' 



(3.16) 



This can be evaluated straightforwardly in dimensional regularization. Adding the contributions from the charged 
pions one has the free nucleon self-energy due to a one-pion loop [58], denoted in the following by T,J, 



uj + iVb 
-u) + i\fh 



327r/2 



(3.17) 



where to = v ■ k = k° and b = — (jj"^ — ie.'^^ In the previous expression we have subtracted the value of the 
one-pion loop nucleon self-energy at w = since we are using the physical nucleon mass. We also need below its 
derivative 



dw 327r2/2 



m^+LO — 3wvo|ilog ^ + TT 



-OJ + I 



Vb 



(3.18) 



The in-medium contribution to the proton self-energy due to the one-7r° loop corresponds to the last line in 
eq. (3.14). Taking also into account the charged pions in the loop we have for the in-medium part of the one-pion 
loop contribution to the proton and neutron self-energies, Sp^^ and ^n,m, respectively. 



U'C 



(27r)4 £2 _ ^2 + 



5{v{k - £)) 



5{v{k - £)) 



0(ej,-|k-l|) + 20(^„-|k-l|) 
e(^„-|k-l|) + 2e(^p-|k-l|) 



(3.19) 



Sij,S^£'^£'' = l2/4 for our choice of v and the second integral on the r.h.s. of eq. (3.14) reads 



(3.20) 



r(i) = o { i^<^(fc°-mep-|k-i|) ^ I d'i i2g(ep~|k-i|) 

J (27r)4 i2-rnl+ie J (27r)3 k^ - P - ml + ie ' 

The step function in the previous integral implies the requirement ^2 > (k _ 1)2 = ]j2 _|_ j2 _ 2|k||l| cos^ . Then, 

k2 -I- l2 _ £2 

cos6'> — ^ =yo ■ (3.21) 



2|k||l| 



It is necessary that yo otherwise cos^ would be larger than 1. This implies that 

|k|-Cp<|l|<|k|+ep- (3.22) 
On the other hand, if |1| > (,p — |k| then j/o ^ ^1- Taking into account these constraints, one has: 
a) |k| > : |1| G [|k| - ^p, |k| + ^ and cos9 G [yo, l] , 

6) > |k| : |1| € [0,^p - |k|] and cosO G [-1, l] ; |1| G [$p - |k|,$p + |k|] and cos0 G [yo, l] . (3.23) 
r(l) 



The same expression of Im results for the cases a) and b). 



7^1) (£p) = -Pl.+ "}^p-y/b arctan^P— P - Vb arctan 



- |k| 



fp+|k| k2-^2_j, 



47r2 



Vb 



Vb 



4|k| 



log 



i^p-Mf + b 



(3.24) 



i^^In all the calculations that follow the square roots and logarithms have the cut along the negative real axis. 
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In terms of Im^ the in-medium part of the one-pion loop contribution to the nucleon self-energy, eq. (3.19), reads 



2 



4/2 

2 

^n,m = ^n,m + ^n,m = (-^^^H?") + ^^m\^p)) ' (3.25) 

where the superscript refers to the pion species in the loop. The full in-medium nucleon self-energy is given by the 
sum of eqs. (3.17) and (3.25). In this way, the proton and neutron self-energies due to the one-pion loop are, in 
that order, 

— ^/ "T ^p,m > 

= + E- „ . (3.26) 
The self-energies for both the proton and neutron can be joined together in E'^, given by 

= i±^E; + ^^E- . (3.27) 

The diagram 5 of fig. 1 originates by dressing the in-medium nucleon propagator, eq. (3.4), with the in-medium 
one-pion loop nucleon self-energy. Its contribution, 115, can then be written as 

n5 = ^Sijk I e*'"" Tr {r" Go(fc) ^^k) Goik)} , (3.28) 

with the convergence factor e^'^"^, 0+, associated with any closed loop made up by a single nucleon line [57]. 
The trace acts in the spin and isospin spaces and gives the result 

n5 = Pijk J (046*'°" {Goik)lE; - Go(fc)^E-) . (3.29) 

Next, we employ the identity 

Goik)l = -^Go{k)i, , (3.30) 
that follows from the r.h.s. of the first line of eq. (3.2). A similar identity also holds at the matrix level 

G,{kf = -^Go{k), (3.31) 

because of the orthogonality of the isospin projectors (1 -|- T3)/2 and (1 — T3)/2. Integrating by parts, as the 
convergence factor allows, we then have 

n» = f « / (g.(.),§ - . ,3.32) 

We perform the integration over k^ making use of the Cauchy-integration theorem. For that we close the integration 
contour along the upper fc'^-complex plane with an infinite semicircle. Because of the convergence factor the 
integration over the infinite semicircle is zero as \mk^ +oo along it. One should then study the positions of the 
poles and cuts in fc° for Go(fc),3 and E^^ in eq. (3.32). First let us note that E/ has only singularities for Imfc° < 0, 
as follows from eq. (3.16). This is also evident for the free part of Go{k)i^, see eq. (3.2). As a result, there is no 
contribution when the integrand in eq. (3.32) involves only free nucleon propagators. The contribution with only 
the density-dependent part both in Go(fc)i3 as well as in the nucleon propagator involved in the loop for E^^ is 
part of the diagram 7 in fig. 1, corresponding to a = 2 contribution. In fig. 3 we depict such an equivalence for 
the diagrams 5 and 7 of fig. 1. An analogous result would hold for the diagrams 6 and 8. The diff'erent Vp = 2 
contributions are evaluated in section 5, so we skip them right now. 
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q 





Figure 3; The equivalence between tiie diagram 5 of fig. f , when only density-dependent parts in the nucleon propagators 
are considered, and the one-pion exchange reduction of the diagram 7 is shown. The second diagram from the left is an 
intermediate step in the continuous transformation of the diagram from the far left to that on the far right. 

Consequently we consider in this section only the contributions where we have simultaneously one free-space 
and one density-dependent part of the nucleon propagators involved in eq. (3.32) and in the calculation of S". Two 
contributions arise. The first one results by employing the density-part for GQ{k)i^ in the r.h.s. of eq. (3.2). The 
integration over /c" is trivial due to the Dirac-dclta function, with the result 



n, 



51 



(27r) 



(3.33) 



feO=£;(k) 



We have introduced the shorter notation 0{£^p — |k|) = 0^ and 0(^„ — |k|) = 0^ . The other contribution involves 




Figure 4: This figure shows that the contribution H^ii, eq. (3.34), involves two free nucleon propagators that follow the 
standard power counting. It is then of 0(jP). 

the free part of Go{k)i^ and can be written as 

ns// = j^sij^ J (2^)4 fco_^(k) + ie Va^^P'™ ~ afc^^^'-j • ^^-^^^ 

However, it is easily seen that this contribution, depicted in fig. 4, is indeed of 0{p^). The reason is because 
there arc two free nucleon propagators of standard counting (those with four-momentum ^ + fc in the figure, 
i — 0(p)), and each of them raises the counting with respect to v given in eq. (2.5) by one power of the small 
scale. As a result, we neglect in the following lis//. The same reasoning is not applicable for lis/ because only 
one nucleon propagator, the one inside the pion-loop self-energy, follows the standard counting. Nonetheless, the 
expression for lis/, eq. (3.33), explicitly shows that it is actually a contribution of ©(p^). This due to the fact that 
dY^jdk'^ = 0{p^), as follows directly from eq. (3.18). We originally counted lis/ as 0{p^) because dJ^J/dk^ was 
taken as 0{p), since T,J = 0{p^) and A;° = O(p^). However, this evaluation of the order of a derivative, based on 
dimensional analysis, represents indeed a lower bound and its actual order might be higher, as it is the case here. 
This mismatch is due to the presence of the variable b, defined after eq. (3.17), in addition to k^. The chiral order 
of the former is fixed by and not by fcg. We also recall that there is another contribution to n5 that results 
by keeping the density-dependent parts both in GQ{k) and S^^, SJ^^^. It will be included in the evaluation of the 
Vp = 2 contributions corresponding to the diagram 7 in fig. 1, section 5. JI5 is an isovector S-wave pion self-energy 
contribution. 
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We consider now the diagrams 6 of fig. 1. The diagram 6a gives 



TTO 

^^6 



9\ 



1 



E{k-q) + i€ 



(3.35) 



where we have omitted the Fermi-sea insertion in the intermediate propagator, fohowing the discussion after 
eq. (3.7). Eg is obtained from Ilg by replacing ^ —q^ and i -H- j in the latter. We now employ eq. (3.30), take 
into account that ct • qa • q = and integrate by parts in In this way eq. (3.35) becomes 



Tja 

'-'-6 



9 9 

4/2 



.9W 



"^'^ Tr [Go(yfc)TV^S^] e^'^"" ^ 



(27r)4 



4/2 J (2^)4 



Tr 



Go(A;)tV-' 



dk° k°-q°- E{Il -q) + ie 



ik°ri_ 



kO-qO - E{k -q)+ie 



(3.36) 



Following an analogous procedure for Hg as the one given below eq. (3.32) for 115, the integration over fc° is 
performed first. As a result, for the previous equation we only take the contributions that simultaneously involve 
one free-space and one density-dependent part of the miclcon propagators in Go{k) and in the loop giving rise to S'^. 
The contribution with only free-space parts vanishes because the integration over fc" along the upper half-plane. 
While that with only density-dependent parts is included in the evaluation of the diagram 8 in fig. 1, section 5. 
On the other hand, applying here the argument in connection with fig. 4, the contribution involving the free-space 
parts of Go(fc) and the density-dependent one of T,^ in eq. (3.35) is 0{p^) because two free nucleon propagators 
(and not just one) are involved with standard counting. Hence, we neglect in the following this other contribution. 
In this way, wc arc left with the contribution, denoted by Ilgj-, that involves T,J and the density-dependent parts 
in Go{k). From cq. (3.36) it is given by 



na 
67 



-I- 



gW [ d^k 
4/2 J (27r)3 

d^k 



4/2 



f _d^ 
J (27r 



Tr 



Tr 



1-T3, 



1 + T3 1 - T3 



ST 



d 



1 



^ dk° k^-qO- E{k -q) + ie 
1 dJlJ 



E(k) - gO - E{li -ci)+ie dk'^ 



feO=B(k) 



feO=E(k) 



(3.37) 



Finally, taking into account the chiral expansion given in eq. (3.8) and adding Ilgj, we have the new quantity Ilg/ 
given by 



n67 = 



f2 qO^'i^J (27r)3^ f 



as, 



dko 



„2 „2 

gA 



k°=E(k) 



n2 r d^k 



(3.38) 



Ilg J, the second 



Hqi is a P-wavc self-energy contribution. However, while the first term on the r.li.s. is isovector 
term is isoscalar, Hgj. It is also the case, see [1], that Hgj is actually one order higher than expected, similarly as 
for S5/. For Hg^ this follows obviously from its explicit expression in eq. (3.38) as SJ = 0{p^). 



In this section we have undertaken the calculation of the diagrams in fig. 1 that can be fully accounted for by 
pion-nucleon dynamics. All the contributions calculated in this sections, Hi to H4, as well as Hs/ and Hg/, are 
linear in density. We have shown that to NLO only the leading contributions. Hi and H2, and the NLO ones H3 
and n4 have to be kept. and Hg/ are finally one order higher. 



4 Nucleon-nucleon interactions 

The inclusion of the nucleon-nucleon interactions for the calculation of the pion self-energy takes place at NLO, 
because Vp = 2 is required at least. As a result, it is necessary to work out the nucleon-nucleon interactions only at 
the lowest chiral order, 0{p'^). These contributions correspond to the diagrams 7-10 in the last two rows of fig. 1. 
First, we discuss these interactions in vacuum and then consider their extension to the nuclear medium. For the 
vacuum case we also discuss nucleon-nucleon scattering up to 0{p). 
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4.1 Free nucleon-nucleon interactions 



The lowest order tree-level amplitudes for nucleon-nucleon scattering, 0{p'^), are given by the one-pion exchange, 
with the lowest order pion-nucleon coupling, and local terms from the quartic nucleon Lagrangian without quark 
masses or derivatives 



^NN = -7;Cs{NN){NN) - -CT{NaN){NaN) . (4.1) 



The fact that these are the leading tree- level contributions is a consequence of our counting eq. (2.5), which 
determines that the lowest order diagrams are those with {di = 0,Vi = l,Wi = 1) and {di = l,Vi = l,Ui = 0). 
The former arises from the contact interaction Lagrangian, eq. (4.1), and the latter corresponds to the lowest 
order one-pion exchange. The tree-level scattering amplitude for iVg^^jj (pi)Afs2,i2 (Pa) ~^ -^83,83 (P3)-^S4, 14 (P4) from 
eq. (4.1) is 

-^ATJV ~ {^3381^3482 ^i3ii^i4i2 ^8382^8481 ^i3i2^i4il) 

Ct {^8381 ' ^8482 ^23*1^2422 ^S3S2 ' ^8481 ^^2322^44*1) 5 (^•^) 

where Sm is a spin label and an isospin one. Obviously, this amplitude only contributes to the nucleon-nucleon 
S- waves. The one-pion exchange tree- level amplitude is 



)('?• q) 8381(0" • q)s4S2 (^4*1 W ■ q')s4si(o'- q') 



NN 



4/2 



S3S2 



+mi- ie q' 



/2 



(4.3) 



with q = p3 — pi and q' = p4 — pi. The corresponding nucleon-nucleon partial waves due to one-pion exchange 
can be calculated using eq. (A. 28). Instead, we first take the one-pion exchange between nucleon-nucleon states 
with definite spin and isospin, so that eq. (A. 28) simplifies to 

N}ne,i,S) = ^l^ E {Oaiai\iSJ){mafai\£SJ)JdpT^;,^{S,I)Yripr, (4-4) 

CTi ,(Tf= — S 

where i and £ are the final and initial orbital angular momentum in the two-nucleon rest frame, respectively. 
Explicit expressions for Nj]{£, £, S) are given in Appendix D of [63]. The sum of the local vertex, eq. (4.2), and the 
onc-pion exchanges, eq. (4.3), is represented diagrammatically in the following by the exchange of a wiggly line, 
fig. 5. 




-I- 



Figurc 5: The exchange of a wiggly line between two nucleons corresponds to the sum of the local and the one-pion exchange 
contributions, eqs. (4.2) and (4.3), respectively. 




Figure 6: Resummation of the two-nucleon reducible diagrams. This is referred in the text as a resummation of the 
right-hand cut or unitarity cut. 

Refs. [3,4] argued that the two-nucleon reducible diagrams should be resummed because they are infrared 
enhanced (by large factors ~ m/|p,|) due to the large nucleon mass. This resummation, depicted in fig. 6, is 
required by our power counting, eq. (2.5), when the latter is applied to the vacuum case as discussed at the 
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Figure 7: Unitarity loop corresponding to tlie function g{A), eq. (4.9). 




Figure 8: Riglit- and Icft-liand cuts of Tji(£,£, S)(p^), for > and < — in order. We liavc also indicated the 
integration contours Ci and Cii used for calculating g{p^) and Njj{£,£, S){p^), eqs. (4.7) and (4.16), respectively. The 
union of both contours Ci U Cii is the one used for Tj"/(£, £, S')(p^), eq. (4.19). In the calculation e — ^ 0"*". 



end of section 2. Notice that every two-nucleon reducible loop in the string is connected by adding O(p^) local 
interactions and the exchange of pionic-lines at the lowest order. As pointed out in the conditions ii) and iii) of 
section 2 the counting does not increase then. Rephrasing the discussion of this section to the present case, the 
nucleon propagators in a two-nucleon reducible loop follow the non-standard counting and each of them is 0{p~'^), 
so that altogether are 0{p~'^). The leading wiggly line exchange is 0{p°). When these two factors are multiplied 
by the 0(p*) contribution from the measure of the loop integrals, associated with the running momenta of the 
wiggly lines, an 0{p^) contribution results. The latter does not increase the chiral order and the scries of diagrams 
in fig. 6 must be resummed.*''^ For this purpose, we follow the techniques of UCHPT [41,42,44] that performs this 
resummation partial wave by partial wave. Many recent nucleon- nucleon scattering analyses using CHPT [5-8] 
follow refs. [3,4] and solve the Lippmann-Schwinger equation in order to accomplish such resummation. UCHPT 
has been applied with great success in meson-meson [42,64,65] and meson-baryon scattering [44,66 70]. The 

#''One could argue that if the nucleoli propagator is taken as 0{p^'^) for the two-nucleon reducible loops, then the measure could be 
taken as 0{p^), counting dp" as 0{p^). If this counting is followed, a suppression by an extra power of p seems to arise. However, this 
factor is multiplied by the large nucleon mass, so that mp finally results, which is then multiplied by local interactions. If the latter 
count as 1/mA, with A ~ m^, the resummation would be required as well within this point of view. We show below that this is the 
case in our approach. 
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master equation for UCHPT is 



Tji % I S)=[I + Njj {£, I S)-g\ Njj {i, l, S) . (4.5) 

This equation, derived in detail in refs. [42,44, 71], results by performing a once-subtracted dispersion relation of 
the inverse of a partial wave amplitude. The function g is defined as follows. Let us denote by p the center-of-mass 
(CM) three-momentum of the nucleon-nucleon system. A nucleon-nucleon partial wave amplitude has two cuts [72] , 
the right hand-cut for oo > > 0, due to unitarity, and the left-hand cut for — oo < p^ < — Tr?4/4, due to the 
crossed channel dynamics. The upper limit for the latter interval is given by the one-pion exchange, as the pion 
is the lightest particle that can be exchanged. These cuts are represented in fig. 8. Because of unitarity, a partial 
wave satisfies in the CM frame that 

ImTj, {£, I S)-} = - , (4.6) 

above the elastic threshold and below the pion production one. The function g in eq. (4.5) only has a right-hand 
cut and its discontinuity along this cut is 2i times the right hand side of eq. (4.6). A once-subtracted dispersion 
relation can be written down given the degree of divergence of eq. (4.6) for p^ oo. The integration contour taken 
is a circle of infinite radius centered at the origin that engulfs the right-hand cut, as shown in fig. 8 by C/. In this 
way 

miA-D) f°° „^ k 

=9iD)-'£{VA-^V\D\) 
.my/A 

^90-1^. (4.7) 

One subtraction has been taken at -D < so that the integral is convergent. Note that the subtraction constant 
g{D) is the value of g{A) &i A = D, in particular, go = g{0)- Since g{p^) = 0{p^), as discussed above, it follows 
that 

go=g{0) = 0{p°) . (4.8) 
The function g{A) corresponds to the divergent integral 

f 1 

The previous integral, depicted in fig. 7, is linearly divergent although it shares the same analytical properties 
as eq. (4.7). In dimensional regularization with _D — > 3 one has, g{A) = —imyfA/ATT. This result is purely 
imaginary above threshold, A> Q, and it corresponds to the imaginary part of eq. (4.7). However, this is just a 
specific characteristic of the regularization method employed, since, as is explicitly shown in eq. (4.7), there is an 
undetermined constant g^. For A — the integral in eq. (4.9) is (infinitely-)negative, so that it is quite natural 
to assume that go < 0. Another more fundamental reason for taking < 0, required by the consistency of the 
approach, is given below. In the following, we regularize any two-nucleon reducible loop in terms of the subtraction 
constant (/o, taking into account eqs. (4.7) and (4.8). The irreducible diagrams with respect to intermediate multi- 
nucleon states will be regularized employing dimensional regularization [73]. This regularization method is shown 
up to NLO in the calculations performed in this work. For explicit calculations of loop integrals apart from g{A) 
within this scheme see Appendix D and the calculation of the energy per nucleon, E / A in section 7. 

Next, we consider how to fix Njj{i, 1, S) in eq. (4.5). This function has only a left-hand cut, due to the exchange 
of pions in the chiral EFT (of course, in a meson-exchange calculation it would include further exchanges of other 
heavier mesons like p, oj, etc). It has no right-hand cut since the latter is fully incorporated in the fimction g{A) 
by construction. As a result, Nji should not be infrared enhanced since the effects of the large nucleon mass, 
associated with the two-nucleon reducible diagrams that give rise to the unitarity cut, are taken into account by 
eq. (4.5). Note that the latter results by integrating over the two-nucleon intermediate states at the level of the 
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inverse of a partial wave, eq. (4.7). In a plain perturbative chiral calculation of a nucleon-nucleon partial wave 
the right-hand cut is not resummcd and the convergence of the perturbative series is spoilt due to the infrared 
enhancement of the two-nuclcon rciduciblc loops. However, since the latter are resummed in eq. (4.5). the idea is to 
match this general equation with a perturbative calculation within CHPT up to the same number of two-nucleon 
reducible loops. The number must be the same to guarantee that Nji is real along the physical region and fulfills 
the requirement of not having right-hand cut. We can make use of the geometric series in powers of g of Tji, 
eq. (4.5), 

rjz(£, I S) = Nji{£, I S) - Nji{e, I S)-g- Nji{£, £, S) 

+ Nji{e, I S)-g- Njiii, I S)-g- Nji{£, 1,8) + ... (4.10) 

where we have used the matrix notation Njj ■ g for the case with coupled channels. Here, g just corresponds to 
the identity matrix times eq. (4.7), because the latter is the same for all partial waves. Together with the previous 
geometric series one also has the standard chiral expansion 

n 

Nji =Y.^n\ (4-11) 

with the chiral order indicated by the superscript. Now, for the determination of the different Nj^\ m < n, the 
matching between eq. (4.10) is performed with a perturbative chiral calculation for which the reducible part of 
every two-nuclcon reducible (or unitarity) loop is counted as 0{p). It is important to stress that this counting is 
applied for calculating Njj not Tjj, for the latter each two-nucleon reducible loop counts as 0{p'^), eq. (2.5). In 
this way, the matching up to a chiral order n automatically comprises at most n two-nucleon unitarity loops. In 
addition, the chiral order of the vertices employed will also make that no spurious imaginary parts are left since 
one is handling in the matching with perturbative unitarity up to order n. 

A few examples will clarify this process of matching and why it makes sense to take as 0{p) the reducible part of 
a two-nucleon reducible loop for calculating Njj within UCHPT. At lowest order, n = 0, there are no two-nucleon 
reducible loops and Nf) {£,l,S) = L^jj{i,i, S), where the latter is the tree-level calculation in CHPT at 0(p°) 
given by the sum of T^^, eq. (4.2), and T^j^, eq. (4.3), projected in the appropriate partial wave. This is the 
wiggly line at the far left of fig. 6. At 0(p), n = 1, the only new contribution is the two-nucleon reducible part of 
the second diagram in fig. 6, denoted by i'yj {£, £, S) for a given partial wave. Writing Nji = Nf^ + TV^ + Oip""), 
and matching eq. (4.10) with the sum of the first two diagrams of fig. 6 one has 

N(f) + _ j^if) . g . Nif) + oip-) = + L^j} + Oip^) , (4.12) 

with the result 

N^'^ = L^j} + Nf^ ■ g ■ (4.13) 
Notice that in the expansion of eq. (4.10) each factor of the kernel Njj{£,£, S) multiplies the loop function g 
with its value on shell. This is why in eq. (4.12) we have —N^j ■ g ■ N^j for one iteration of g, which is then 
subtracted from the function L^jj in eq. (4.13). This equation shows explicitly that the simultaneous expansion in 
chiral powers and number of loops for fixing N^jf^ implies that UCHPT really takes as 0{p) the difference between 
a full calculation of one two-nucleon reducible loop and the result obtained by factorizing the vertices on-shell, 
eq. (4.10). Ultimately this relies on the fact that the difference has no right-hand cut, which is the one associated 
with the infrared enhanced two-nuclcon reducible loops, and it has only left-hand cut. The latter is incorporated 
perturbatively in the interaction kernel Nji{£,£, S), which is improved order by order. This is the reason why we 
have treated the expansion in two-nucleon reducible loops on the same foot as the chiral expansion. This procedure 
is iterated up to any desired order. E.g. at O(p^) new contributions would arise that require the calculation of 
the irreducible part of the box diagram in fig. 6 and the reducible parts of the last diagram of fig. 6 with the 
wiggly line exchange iterated twice [74] . In addition, there are also local interaction terms from the quartic nucleon 
Lagrangian and two-nucleon irreducible pion loops [7, 8, 25, 73]. If we denote^ all these new contributions projected 

(2) 

(into the corresponding partial wave by L jj{£,£, S), the following equation results 

Nf^ = L + AT « . , . AT + AT • , • TV « - Nf^ ■ g ■ N^^ ■ g ■ Nf^ . (4.14) 
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The Nji calculated up to some given order in the expansion eq. (4.11) is then substituted in eq. (4.5). On the 
other hand, one can match formally eq. (4.5) with a perturbative chiral calculation of Tji{i,l,S) for any value 
of the S-wave nucleon-nucleon scattering lengths because the latter enter parametrically in the calculation. This 
procedure gives rise to values of the low-energy constants Cs and C't that are consistent with their ascribed 0(p") 
scaling, see eq. (4.28) below. It is worth pointing out that eq. (4.5) is algebraic, so that the numerical burden for 
in-medium calculations is reduced tremendously. 




Figure 9: Box diagram, L^Jj, originating from the first iteration of a wiggly line. It consists of the diagrams shown on the 
right-hand side of the figure with zero, one or two local vertices and/or one-pion exchange ampUtudes. 

The dependence on the parameter go takes places because of the infrared enhanced two-nucleon reducible loops. 
This has made necessary to rcsiim the right-hand cut, which requires the presence of one subtraction constant, go, 
eq. (4.7). Indeed, for a fixed chiral order, according to the application of eq. (2.5) to nucleon-nucleon scattering in 
vacuum, the dependence on go becomes smaller as higher powers of g are considered for calculating Njj, eq. (4.11). 
To show this, we need to take advantage of the analytical properties of iVj/, cq. (4.5). As discussed above, this 
quantity has only the left-hand cut, shown in fig. 8. For the following discussion we take the case of one uncoupled 
channel to simplify the writing. Its generalization to coupled channels is straightforward employing a matrix 
notation. The imaginary part of Njj along the left-hand cut is given by, 

ImNji = \^^lmTjj = |1 + gNji\'lraTji , |pp < . (4.15) 

Note that g is real along the left-hand cut. We employ this result to write down a once-subtracted dispersion 
relation for Njj. The integration contour is shown in fig. 8 as C// and consists of a circle of infinite radius centered 
at the origin that engulfs the left-hand cut. 

We have taken one subtraction in the dispersion relation because the one-pion exchange amplitude, eq. (4.3), tends 
to a constant for — ^ oo. Then, we have for Tjj, eq. (4.5), 



Tji{A) = 



A-D lmTji\l+gNjj\^ 



Nji{D) + / dk 



{k^ - A - ie){k^ - D) 




(4.17) 



In order to solve eq. (4.16) one needs ImTji as input along the left-hand cut. CHPT could be used, since this 
imaginary part is due to multi-pion exchanges. As a result one could afford its calculation perturbatively because 
the infrared enhancements associated with the right-hand cut are absent in the discontinuity along the left-hand 
cut. The reason is because this discontinuity, according to Cutkosky's theorem [75,76], implies to put on-shell 
pionic lines so that within loops the pion poles are picked up making that the energy along nucleon propagators 
now is of 0{p), instead of a nucleon kinetic energy. In this way, the order of the diagram rises compared to that 
of the reducible parts and it becomes a perturbation. E.g., let us take as illustration the last diagram on the r.h.s. 
of fig. 9, corresponding to the twice iterated one-pion exchange. Its reducible part is infrared enhanced, which 
has been calculated by us in the presence of the nuclear mcdiiim. in agreement with ref. [59] when reduced to the 
vacuum case. However, its discontinuity across the left-hand cut arises by putting on-shell the two intermediate 
pion lines. Its leading contribution to ImT in a 1/m expansion in the f-channel CM frame is given by 

lrnT=-M ^'/'-jJ^''\ t>4ml, (4.18) 
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with t — — 2p^(l — coaO), cos9 G [—1,1] and A/" is a numerical factor due to the spin algebra. No factor m 
appears in the numerator and it follows the standard chiral counting. In this way, the leading contribution to 
ImTji along the left-hand cut is given by the one-pion exchange. The latter can then be inserted in eq. (4.16), 
once projected in a given partial wave. The solution of this equation would correspond to the leading result for 
Nji in the chiral expansion of eq. (2.5), without involving the expansion in the number of two-nucleon reducible 
loops. This interesting exercise will be left for future consideration. 

Pion exchange amplitudes are treated perturbatively in the Kaplan-Savage-Wise (KSW) power counting [12, 13]. 
This is done for any energy region and, in particular, along both the right- and left-hand cuts. On the other hand, 
the dispersive treatment offered here only needs as input the discontinuity (imaginary part) of a nucleon-nucleon 
partial wave along the left-hand cut, see eq. (4.17). This discontinuity arises due to pion exchanges which, as 
discussed in the previous paragraph, could be calculated perturbatively in CHPT. Differences with respect to 
KSW arise due to the resummation of the right-hand-cut in eq. (4.17), including both local and pion-exchange 
contributions. This would correspond to higher orders in KSW power counting [13]. Notice also that while KSW is 
a strict perturbation theory calculation in quantum field theory (QFT) ours merges inputs from perturbative QFT 
and S-matrix theory, see e.g. ref. [77] for a pedagogical account of first applications of the similar N/D method to 
nucleon-nucleon scattering. 

Two subtraction constants appear in eq. (4.17), Njj{D) from eq. (4.16) and go from the function g{A), eq. (4.7). 
We are going to show that they are not independent, however. The two constants have appeared due to the splitting 
between the functions Njj and g when expressing Tj/ = NJ^ + g, eq. (4.5). This is analogous to the standard fact 
that in any renormalization scheme there is an exchange of contributions between local parts in loops and local 
counterterms. In order to proceed with the demonstration that the resulting T//, eq. (4.17), does not depend on 
the subtraction constant g{D), let us write directly a dispersion relation for TJ"/ taking the contour C/ U C// in 
fig. 8 

m{A-D) k {A-D) /""''/'^^a ImTji l\Tji\'' 



J JI {A) - 1 jj [U) I * {k^-A- ie){k^ -D) TT y_ 



(fc2-A-e)(A;2-£>) 

A-D S-^ fjiiP) 



OO 

(4.19) 

fe2=0 



{I - 1)! d(fc2)^-i (A;2 - A)(fc2 - D) 



where fjiip"^) — \p\'^^Tjj^{p'^) . The last term in the previous equation gives contribution for I > 1 and arises 
due to the behaviour at threshold of a partial wave, vanishing as [pl^^. Two-body unitarity is assumed all the way 
along the right-hand cut in the first integral. This is not essential for the discussion that follows and we could have 
written directly ImTJj^ along the right-hand cut, as done for the left-hand one."^^ As discussed above the input 
for solving Nji in eq. (4.16) is ImTjj along the left-hand cut. This can also be shown explicitly from eq. (4.19) 
by writing l/jTj/j^ = \Njj + gj^, as follows from eq. (4.5). Subtracting g from TJ"/ we arrive to the following 
equation for Njj, 



(4.21) 



{k^-A-e){k^-D) 

A-D fjj{k^) 



{£ - 1)! d(A;2)^-i (P - A){k^ - D) 



fe2=0 



In the following we omit the last term in the previous equation for simplicity, since it does not depend on g{D). The 
reader could include it straightforwardly if desired. If eq. (4.21) is solved by iteration, it is straightforward to show 
that Tji does not depend on go at any order in the iteration. The zeroth iterated solution is NJj.q = Ty/{D)—g{D), 



#4 In eq. (4.19) we could use different subtraction points for the two integrals, e.g. B and D, respectively. One then has 
Tji (A) - Tjj (D) + dk .^-^—-^-^—^ J^^ dk (^i^2_A-ie){k^-B) 



A-D f-"''/\^2 lmTjj/\Tjj\^ A-D d'-^ fji{k^) 



{k^-A- ie){k^ -D) {e - 1)! d(fe2)«-i (ifc2 _ A)(fe2 _ D) 



(4.20) 
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which yields TjIq{D) = TJ^iD) - g{D) + g{A). Obviously, the sum —g{D) + g{A) is independent of g{D). For 
the first iterated solution one has 

Notice that only the combination —g{D) + g{k'^) appears in the integral, which is independent of g{D). However, 
^jhi depends explicitly on g{D) due to the term before the integral. Nevertheless, given that TJj.^{A) = Njl.^{A) + 
g{A), the first g{D) on the r.h.s. of eq. (4.22) is accompanied again with g{A) so that no dependence on g{D) 
is left. This process can be straightforwardly generalized to any order. For the j*^ iteration the combination 
T~^{D) — g{D) that appears in Nj^.j{A) before the integral is added to g{A) for calculating Tji-j{A), so that no 
dependence on g{D) arises from this fact. In addition, under the integration sign we have repeatedly j times the 
same term Ty^{D) — g{D) + g(k'^), which does not depend on g{D). 

From the previous discussion, one concludes quite confidently that no g(D)-depcndence is left because this was 
the case for Tjj{A) evaluated at any order in the iterative solution of Njj{A), eq. (4.21). In this way, it is clear that 
one could interpret the constant g{D), eq. (4.7), and the subtraction point D in close analogy with renormalization 
theory. The latter corresponds to the "renormalization scale" and the former fixes the "renormalization scheme" . 
For a given g{D) then Njj{D) is fixed so as to reproduce Tjj{D) at the point |p^| = D. The dependence on g{D) 
is then transmuted into the experimental input Tji{D). The final result should be independent of g{D), which in 
turn, by taking the derivative of TJ"/, eq. (4.17), with respect to this parameter implies the equation 

dNjijP) 2 2{A - D) ImTji Re [{gdNjj/dg{D) + Nji){l + N*,g)\ 

-^m- = ^ i-oo * [k^-A-ie)[k^-D) • ^'-^'^ 

This discussion also shows that one always has the freedom to take go to be the same for all the partial waves, as 
we have done.'^^ 

For higher partial waves it is convenient to derive the dispersion relation for iVj// |p|^^ instead of eq. (4.17). 
In this way, the low energy behaviour of a partial wave as |p|^^ for |p| — )• is ensured, independently of the 
approximation for ImTj/ [80]. The resulting expression is 

Note also that for ^ > 1 no subtraction is needed if ImT const. (mod log) for fc^ oo, as in the one-pion 
exchange. Then, one could also rewrite the previous equation for ^ > 1 as 

The degree of divergence of ImTj/ for |p| — >■ oo increases by including higher order loop contributions, see e.g. 
eq. (4.18). As a result, more subtractions should be taken and the resulting subtraction constants could be related 
with higher order chiral counterterms. 

We have proposed to consider g as 0{p) in order to fix Njj. Indeed, g is suppressed along the left-hand cut, 
vanishing in the low momentum region of the dispersive integral of eq. (4.16), which dominates its final value for 
low energy nucleon-nucleon scattering. On the physical Riemann sheet |k| = +iK, with /t = \/ —k"^ > 0, and since go 
is negative and of natural size — mm7r/47r, it tends to cancel with — im|k|/47r = mK/iw > and becomes zero for 



#''Thcrc is an infinity of solutions of eq. (4.19) differing between each other in the number of zeros of Tjj. Each of these zeros is a pole 
of TJj so that it brings altogether as free parameters the position of the pole and its residue. They are the so-called Castillejo-Dalitz- 
Dyson (ODD) poles [78]. The ODD poles are typically associated with resonances [78,79]. Notice that a pole in TJ^ typically maJjes 
its real part to vanish if the remnant is a smooth function of energy around the pole. In low energy S-wave meson-meson scattering 
the Adler zeros correspond to ODD poles [42]. However, for nucleon-nucleon scattering there is no evidence for a low energy zero in 
the partial waves (apart from the trivial one at threshold for i > 1.) In the pionless EFT for nucleon-nucleon interactions the third 
integration on the r.h.s. of eq. (4.19) is absent. The infinity tower of chiral counterterms in this EFT can be accounted for by adding 
CDD poles, see ref. [42] where this is shown explicitly for a similar problem. 
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K = —Ango/m ~ rnvr- This is an important reason for having taken go < above. Proceeding along these lines, so 
that g is treated as relatively small along the left-hand cut, eq. (4.16) would simplify at leading order. On the one 
hand, |1 + gNji\'^ is replaced by 1 and, on the other, ImT// is given by the onc-pion exchange. Hence, one obtains 
for iV(°; in S-wave the sum of a constant plus one-pion exchange (resulting from the dispersive integral), precisely 
the content of the wiggly lines, fig. 5. For ^ > 1 let us take directly eq. (4.25). In this way, when neglecting gNjj, 
the dispersive integral just gives rise to the onc-pion exchange, as was the case for our previously calculated N^jf . 
One could continue further in this way, and solve eq. (4.16) in a power series expansion of g along the left-hand cut 
at each chiral order in the calculation of ImTjj. The truncation of such expansion leaves a residual go dependence. 
We have followed the same point of view in order to determine Nji through the matching process discussed above. 
Indeed, alternatively to performing the geometric series expansion of eq. (4.10), we could consider directly the 
inverse of Tjj, similarly as done in order to obtain eq. (4.15). Then, it follows from eq. (4.5) that 



T* N* 



+ 9 

2 IAr_j2 



Nji = Tji\l + gNjil' - \Nji\'g* . (4.26) 

The first method discussed above for determining Njj is the perturbative solution of eq. (4.26) in a chiral series 
of powers of g. The solutions of eqs. (4.26) and (4.16) employing the perturbative method are equivalent because 
Nji from eq. (4.26) has only a left-hand cut, being its imaginary part along this cut the same as eq. (4.15), and it 
is analytical, so that it satisfies the perturbative version in power of g of the dispersion relation eq. (4.16).*^ The 
following remark is in order. The perturbative solution in the chiral expansion of powers of g of eq. (4.26) has the 
advantages over solving eq. (4.16) that it is algebraic and the chiral countcrtcrms in Tji arc taken into account in the 
solution Njj in a straightforward manner. It is also very versatile, so that it can be extrapolated straightforwardly 
to correct by initial and final state interactions and to the nuclear medium. Notice that eq. (4.26) can only be 
solved pcrturbatively since the input Tji is calculated in CHPT and only fulfills unitarity perturbatively. However, 
the exact solution of the integral equation eq. (4.16) has the advantage of not requiring the expansion in powers 
of g but just the chiral series on ImTjj along the left-hand cut, and the latter expansion rests in a sound basis as 
discussed above. 

We now concentrate on fixing the constants Cs and Ct from the local quartic nucleon Lagrangian, eq. (4.1). 
These constants and go, eq. (4.8), are the only free parameters that enter in the evaluation of the nucleon-nucleon 
scattering amplitudes from eq. (4.5) up to 0{p). We first discuss the LO result and then the NLO one. Cs and 
Ct are fixed by considering the S-wave nucleon-nucleon scattering lengths a* and for the triplet and singlet 
channels, respectively. At 0{p^) we have at threshold 

ror(0,0,0)= -^^^-3^^) 



l-goiCs-SCr) ' 

The triplet S-wave is clastic at this energy, without mixing with the ^Di partial wave, because of the vanishing of 
the three-momentum. The resulting expressions for the scattering lengths from eq. (4.27) imply that 

m Wngo/m + 3/as + 1/at 
IGtt {go + m/{4Tras)){go +m/{'i'!Tat)) 

Ct = — ^/'^s-l/at _ .4_28) 

167r {go + m/{ATras)){gQ + m/(47rat)) 

One of the benchmark characteristics of nucleon-nucleon scattering are the large absolute values of the S-wave 
scattering lengths = — 23.758±0.04fm and at = 5.424±0.004fm, sothat m^r '> \'i-/as\, l/^t- Given the expression 



i^^We have shown this equivalence explicitly for N^) . It is also straightforward to show it for N^) . 
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for the imaginary part of g{A) above threshold in eq. (4.7) one can estimate that go ^ —mrriTr / 4tt ^ —0.54 m^,'^^ 
which is then much larger in absolute value than m/(47r|as|) and m/(47ra(), although there is a difference because 
at is smaller by around a factor 4 than \as\- As a result, it follows from eq. (4.28) that \Cs\ ~ l/|ffo| ^ \Ct\ = 
0{m/16TratgQ). In this way, the low-energy constants Cs and C't do not diverge for Us, a* — > c» and after iteration 
it is still consistent to treat 'Cj^^y: ^q. (4.1), as 0{p^). Notice as well that the one loop iteration of the contact 
terms compared in absolute value with the tree level goes like —m\p\{Cs — (45'— l)CT)/47r, taking into account the 
expression for g{A) given in eq. (4.7). The three-momentum is divided by the scale ~ —An/mCs ~ ttIt^, considering 
the just given estimates for Cs ^ l/go and go ^ — TOm,r/47r. This justifies to iterate these diagrams for |p| = 0{p) 
as discussed above. For the case of the once-iterated pion exchange one would have the factor m\Tp\g\/ l&ir f"^ as 
compared with the tree level one-pion exchange. Then |p| is divided by the scale l%-Kf^/mg\ ~ 2m.,^ = 0{p), and 
the one-pion exchange should be as well iterated together with the lowest order contact terms. The issue of iterating 
potential pions is analyzed in detail in ref. [13], in order to understand the failure of the KSW power counting in 
some triplet channels, particularly, for the ^Si—^Di and ^Po,2 channels. The authors of ref. [13] conclude that for 
some spin triplet channels the summation of potential pion diagrams is necessary to reproduce observables, while 
for the singlet channels this iteration does not seem to be a significant improvement over treating pion exchanges 
perturbatively. 




go (jnl) 

Figure 10: (Color online.) Values for li (red solid line), So) (magenta dashed line) and £2(^Si) (cyan dot-dashed line) 
as a function of go- ^2 is expressed in units of . 

Only local terms and one-pion exchange contributions enter in the calculation of N^j {£,£, S). This is rather 
simplistic in order to describe properly the nucleon-nucleon interactions as a function of energy soon above threshold. 
Let us now consider eq. (4.5) with Njj up to 0{p). At this order, ImTji along the left-hand cut is stiU given by the 
one-pion exchange, so that the exact solution of eq. (4.16) for Njj would be the same. The differences observed in 
the results at 0{p^) and 0{p) are then due to keep a one more factor g in the perturbative solution of eq. (4.16). 
This discussion shows clearly the mixed nature of the chiral expansion in powers of g for obtaining Njj. 

We employ the ^5*0 and ^Si scattering lengths for evaluating Cs and Ct at 0{p). We denote by a any of these 
scattering lengths and apply eq. (4.5) at threshold. We obtain 



1 ImTji 



kReTji 



m 



fc-i-O 



Nji 



fc-i-O 



(4.29) 



Taking eq. (4.13) at threshold we rewrite Nf'j = —C and express L 



(1) 
// 



-C'^go + C£i +£2 because the box diagram 
l'^j], fig. 9, consists of four contributions with two, one and zero local vertices. The first contribution is given by 
—C^go, the second by C£i and the last one by £2, respectively. The coefficients £1 and £2 are given in terms of go 
and the known parameters m, gA and m-n-. £1 is the same for the partial waves ^So and ■^S'l while £2 is different. 



*7As explicitly shown in the second line of eq. (4.7) one can trade between the subtraction constant and 
the subtraction point. In a natural way, both should be taken of similar size for estimations. 



/~A/4:TT just by changing 
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The values of £i and £2 as a function of go are shown in fig. 10. Substituting these expressions in eq. (4.29) 



C-^, (4.30) 



with = + go) the O(p0) result. 
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Figure 11: (Color online.) ^So, '^Si, "^Di phase shifts and the mixing angle ei as a function of |p|. The (red) solid lines 
correspond to the Nijmegen data [81,82]. For the rest of the lines three values of go = —1/4, —1/3 and —1/2 are 
employed. For LO these lines are the (green) dashed, (cyan) dot-dashed and (magenta) dotted lines, respectively. While 
to NLO these are the (black) double-dotted, (orange) double-dot-dashed and (blue) short-dashed lines, in that order. The 
squared points in the panel for ei corresponds to a NLO calculation with go = —0.1 m^. 

In figs. 11, 12 and 13 we show the LO and NLO results for the nucleon-nucleon scattering data (phase shifts 
and mixing angles) up to |p| = 300 MeV making use of eq. (4.5). Since Cs at LO is close to l/go, as explained 
above, we show the results for the values go = —1/4 m^, —1/3 and —1/2 because its inverses are —4 rn~^, 
—3 m^^ and —2 m^^, respectively. In this way, the resulting Cs at LO is of order 1 m^^, a natural size. E.g. 
employing the estimation for go — —0.54 given below eq. (4.28), one would obtain l/go ^ —1.8 m~^. On the 

other hand, let us recall that negative values for go, and not far from —0.5 mj, are the required ones in order to 
optimize the perturbative solution of eq. (4.16). For the LO results the lines are the dashed, dot-dashed and dotted 
lines, corresponding to go — —1/4, —1/3 and —1/2 m^, respectively. While to NLO these are the double-dotted, 
double-dot-dashed and short-dashed lines, in the same order. For |p| ~ 360 MeV the pion production threshold 
opens and it does not make sense to compare with data above this point. An 0{p^) calculation, which includes 
important new physical mechanisms, as non-reducible two-pion exchanges between others, as indicated above before 
eq. (4.14), is presumably needed to improve the agreement with data [5,8]. E.g., it is well known that for the ^5*0 
partial wave an 0{p^) chiral counterterm, in the standard chiral counting, is required in order to reproduce its 

#8 At 0{p°) in the KSW counting [12]. 
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Figure 12: (Color online.) ^Pi, ^Po, ^Pi, ^^2, phase shifts and the mixing angle £2 as a function of |p|. For notation, 
see fig. 11. 



relatively large effective range so that the agreement with data improves. This can be understood by considering 
the effective range expansion. For the ^ 5*0 partial wave l/a^ is extremely small so that the contribution from the 
effective range rop^/2 rapidly overcomes — l/a^ (the leading order contribution). Then, this problem is not so 
much related to the fact of having too large higher order corrections but more it arises because the leading order 
is anomalously small. The largest differences in absolute values between the LO and NLO results are observed in 
the ^Si-^Di and ^Po partial waves. These partial waves, as discussed in depth in ref. [13], have large non-analytic 
corrections from two potential pion exchange. For the ^Pi, and ^D^^ waves the difference in absolute terms 
is small, a few degrees, although relatively it can be large typically for |p| > 150 MeV. For higher partial waves 
these differences are typically much smaller since the iteration of one-pion exchange becomes smaller [59]. Our 
0{p) results are of comparable quality to those obtained at LO within the Weinberg's counting approach [8]. The 
^Pq phase shifts are also not well reproduced at this order in ref. [8]. Both approaches share the same input for 
ImTj/ along the left-hand cut, and at 0{p) we have already considered the iteration of one g factor in determining 
Nji, as discussed above. The main differences between our results and ref. [8] at LO concern ti and the phase 
shifts for ^Pi and ^D^. For the latter our results are closer to experiment while for the two former observables 
the LO calculation of ref. [8] is closer to data. It is known that one-pion exchange has a too large tensor force 
which is reduced by higher order counterterms. In the meson exchange picture this cancellation at short distances 
of the one-pion exchange tensor force is produced by the exchange of p-mesons [83]. The mixing '^Si-^Di and the 
partial wave ^Pq have large attractive matrix elements of the one-pion exchange tensor operator, as stressed in 
refs. [15, 21]. These are the partial waves that depart more from data in absolute terms. The '^Pq phase shifts were 
reproduced accurately in ref. [15] at LO for low energies. In this reference, a counterterm was promoted to LO in 
all the partial waves with attractive tensor interactions, and in particular to the ^Pq channel. Such free parameter 
is needed to fit the '^Pq data [8, 15]. The results of ref. [15] are then cut-off independent for high enough values of 
the employed cut-off. 
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Figure 13: (Color online.) ^-D2, ^^32, ^fa, ^Ga phase shifts and the mixing angle es as a function of |p|. For notation, see 
fig. 11. 



As a result of the perturbative approach actually followed in this paper for determining N ji by solving eq. (4.26) 
in an expansion in the number of two-nucleon reducible loops, a residual dependence on go is left in the solution 
due to higher orders in this expansion (and not from the pure chiral one, eq. (2.5)). As more orders are included 
the exact solution of Tj/, obtained by solving eq. (4.16), is better approached and any dependence on go should 
tend to vanish. From here one could also infer that contributions with one-pion exchange twice iterated in iVj/ 
are expected to be significant at least in those observables with a clear go dependence in figs. 11-13. It is also 
worth noticing that the dependence on go in figs. 11-13 at LO/NLO is much smaller for the P- and higher partial 
waves than for the S-waves. This should be expected because Nji for > 1 vanishes at threshold as |pp^ so 
that both g and iVj/ are small in the low energy part of the left-hand cut. In this way the perturbative solution 
of eq. (4.16) should typically converge faster for higher Conversely, the convergence of the S-waves should be 
slower, something that it is clear for the Di coupled channels from fig. 11. Particularly noticeable is the 

dependence on go of ei, a fact that is in agreement with the results of Fleming, Mehen and Stewart [13]. The 
squared points in the panel for ei in fig. 11 are obtained with go = —0.1 mj. They agree closely with data [81], 
though such good agreement seems to be accidental. 



4.2 Nucleon-nucleon scattering in the nuclear medium 

When calculating a loop function in the nuclear medium we typically use the notation Ly , where i indicates the 
number of two-nucleon states in the diagram (0 or 1) and j the number of pion exchanges (0, 1 or 2). In addition, 
we also use Lijj, Lij^m and Lij^d, with the subscripts /, m and d indicating zero, one or two Fermi-sea insertions 
from the nucleon propagators in the medium, respectively. In this way, the function g — Liqj and its in-medium 
counterpart is iio, that is calculated in the Appendix C. 

The evaluation of the nucleon-nucleon scattering amplitudes in the nuclear medium at lowest order can be 
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easily obtained from our previous result in vacuum since the only modification without increasing the chiral order 
corresponds to use the full in-medium nucleon propagators. This is directly accomplished by replacing g{A) by Lio 
in eq. (4.5). At any order for nucleon- nucleon scattering in the nuclear medium, we use eq. (4.5) but now with the 
function g substituted by Lio so that 

T^jii, i, S)=[I + N%{1, 1 S) ■ L%] . Nyj{£, £, S) . (4.31) 

The same process as previously discussed is followed to fix Njj. Note that any other in- medium contribution 
requires Vp = 1, which then increases the order at least by one more unit, cf. eq. (2.4). This new in-medium 
generalized vertex must be associated with the nucleon-nucleon scattering diagrams of leading order. The modifi- 
cation of the meson propagators (both heavy and pionic ones) by the inclusion of an in-medium generalized vertex 
increases the chiral order by two units. However, the modification of the enhanced nucleon propagators with one 
in- medium generalized vertex only increases the order by one unit and these contributions must be kept at NLO. 
It goes beyond the scope of this article to offer a complete study of the in-medium pion self-energy at N^LO where 
the full NLO in-medium nucleon-nucleon interactions are needed. What we do here for illustration is to change the 
free nucleon propagators by the in-medium ones in the calculation of the box diagram L^j] that enter in fixing A^j^ , 
eq. (4.13), with g also replaced by Lio. In eq. (4.31) we have included the superscript ■is, which corresponds to 
the third component of the total isospin of the two nucleons involved in the scattering process, both in the partial 
wave Tyj{£,£, S) and in L^q, as the Fermi momentum of the neutrons and protons are different for asymmetric 
nuclear matter. The function L'-^q conserves total isospin /, because it is symmetric under the exchange of the two 
nucleons, though it depends on the charge (or third component of the total isospin) of the intermediate state. This 
is a general rule, all the is = operators are symmetric under the exchange p O n, so that they do not mix isospin 
representations with different exchange symmetry properties. 

In this section we have determined the vacuum nucleon-nucleon scattering at LO and NLO following the novel 
counting of eq. (2.5) [1]. For nuclear matter the LO nucleon-nucleon scattering amplitudes have been also obtained. 
The infrared enhancement of the two- nucleon reducible loops have made it necessary to resum the right-hand cut. 
This is accomplished by a once-subtracted dispersion relation of the inverse of a partial wave giving rise to eq. (4.5). 
The important function g{A), eq. (4.7), which is defined in terms of a subtraction constant, g{D) or go, is intro- 
duced. It has been argued that the subtraction constant is 0{p°), because by changing the subtraction point B 
the subtraction constant is modified reshuffling the form of the function g{A), which is invariant. The process 
for determining the interaction kernel Njj, eq. (4.5), has been also discussed in detail. It was obtained that the 
subtraction point D acts as a "renormalization scale" where an experimental point is reproduced. The subtraction 
constant g{D) just fixes the "renormalization scheme" and the exact results should not depend on it. A natural 
value for go ^ —mniT^/A'K was argued to be adequate for obtaining Nji as a perturbative solution of eq. (4.16) in 
order to suppress the effects of the iterative factor \l + gNji\'^ in the equation. The couplings Cg and Ct from the 
local nucleon-nucleon Lagrangian, eq. (4.1), have been fixed in terms of go at LO and NLO reproducing the S-wave 
nucleon-nucleon scattering lengths. These couplings keep their estimated size of O(p^) after the iteration, despite 
the well known fact that the nucleon-nucleon scattering lengths are much larger than l/niTr- The resulting phase 
shifts and mixing angles at LO and NLO are depicted in figs. 11, 12 and 13. It is argued that higher orders should 
be included in order to improve the rcprodiiction of data. Particularly, a N^LO analysis shoiild be pursued since it 
would include the important two-pion irreducible exchange and new counterterms, in particular the one necessary 
to reproduce the effective range for the ^So partial wave [11]. This is left as a future task since our present main 
aim is to work the results up to NLO and settle the formalism in detail. 

5 Contributions from the nucleon self-energy due to nuclear interac- 
tions 

In this section we consider those diagrams in fig. 1 that include the nucleon-nucleon contributions to the nucleon 
self-energy in the medium, diagrams 7 and 8. In turn, for each of these figures the one on the top corresponds 
to the direct nucleon-nucleon interactions, while the exchange part gives rise to the diagram on the bottom (that 
includes the part of the diagrams 5 and 6 with all nucleon propagators corresponding to Fermi-sea insertions.) 
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First, let us consider the evaluation of the diagrams 7 in fig. 1, denoted by Ily. It is given by 

ny = ^eijk J ^e"'"" Tr {T>'Go{h)J:NNGo{h)} , (5.1) 
where ^ ^ 

SjVAT = ^ Sp^jvJV H ^ ^n,NN , (5.2) 

with Sp^ATjv and S„,jvAr the proton and neutron self-energies due to the nucleon-nucleon interactions, in order. 
Performing the trace in isospin, 

nr = ^£.,3 E / (Go(fci)^Sp,iViv - Go(fci)^S„,^iv) • (5.3) 



Here ai corresponds to the spin of the incident nucleon. Taking into account the identity eq. (3.30) we can integrate 
by parts eq. (5.3) with the result 

TT f ^^^'^ ik% fry II \ 9^p,NN ^ ^ 5E„ Arjv\ 




Figure 14: In-medium nucleon self-energy due to the nucleon-nucleon interactions with the Fermi-seas. 
The nucleon self-energy due to the nucleon-nucleon interactions, represented in fig. 14, is given by the expression 



''ai,NN — ~^ 

o:2if 2 



— -^e**'°''Go(fc2)a2^JViv(fclCriai,fc2Cr2Q;2|fci(Tiai, A;2(T2CI!2) (5.5) 



where Tat^v is the two-nucleon scattering operator between the nucleon states characterized by the four-momentum 
ki , spin ai and third component of isospin ai . We also use the variables 

a=^{ki+k2) , p=^{ki-k2) , (5.6) 

and 

A = 2ma° - , (5.7) 

with a the three-momentum made up from a% i = 1,2, 3. We introduce the shorter notation 

^aiS (P' =TNN{kiaiai,k2(72a2\kiaiai,k2a2a2) , (5.8) 

that is more convenient for forward scattering than the notation followed in Appendix A. For on-shell scattering 
A = p^. Eq. (5.4), after using eq. (5.5), becomes 

^7--^2^e,,■3^^^y ^"e \^Go{k,),Go{k2)p 

- Go (fci )„Go (fc2)„ ^^"%fc'o' j • (5-9) 
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In order to obtain this result we have used that 



E 



dki dk2 ik% ik° 
(27r)4(27r)^ 



d 

Go{ki)pGo{k2)n-g^TNN{ki(Tip, k2<J2n\ki(Jip, k2<J2n) 
d 

- Go{ki)nGo{k2)p-g^TNN{kiain, k2(J2p\k\(J\n, k2<J2p) 



(5.10) 



which fohows for two reasons. First, let us notice that because of Fermi-Dirac statistics 

TNN{kicrip, k2cr2n\kiaip, k2cr2n) = TNN{k2cr2n, kiaip\k2(T2n, kiaip) . (5-11) 

Second, at LO the amplitude T/v^r, as commented above, is given by the iteration of the wiggly line in fig. 6. The 
latter docs neither depend on fc^ nor on k^, see eqs. (4.2) and (4.3). Since V^q depends on k^ and fc^ only through 
their sum, A:^ + fc^, then T^n at LO only depends on them in the same way and dT^N/dk^ = dTMN/dk2 holds. 
Taking these two facts into account, as ki and at are dummy variables, eq. (5.10) is obtained. 

It is convenient to give the nucleon-nucleon scattering amplitude as an expansion in partial waves, eq. (A. 8). 
The partial wave decomposition of the nucleon-nucleon amplitudes is derived in detail in Appendix A. A nucleon- 
nucleon partial wave is denoted by Tj^j {£',£, S), where J = £ + S is the total angular momentum, / is the total 
isospin, is = ai+a2, £' and i are the final and initial orbital angular momenta, respectively, and S is the total spin. 
The partial wave is a function of , and A for our previously calculated nucleon-nucleon amplitudes. Since for 
our present case, cq. (5.9), A ^ p"^ an analytical extrapolation in A of T'jj{£', £, S) is necessary. While eq. (A. 8) is 
given in the CM of the two nucleons involved in the scattering process, eqs. (5.1) and (5.5) are given in the nuclear 
matter rest-frame. This implies that one must take into account the boost from the former frame to the latter in 
order to use eq. (A. 8). However, as is shown in Appendix C of ref. [63], the angle of the associated Wigner rotation 
is suppressed and it is C((p/m)^). Then, the leading and next-to-leading nucleon-nucleon scattering amplitudes 
can be used as Lorentz invariants, similarly as for the meson-meson ones, and eq. (A. 8) can be directly used in 
eq. (5.5). Let us recall that our calculation of the pion self-energy in nuclear matter is up to NLO, 0{p^), and these 
relativistic corrections are of 0{p''). From eqs. (5.3) and (5.5) one has to sum over the spins ai and a2- The fact 
that both the initial and final nucleon-nucleon states are the same implies a great simplification in the equations. 
First, if we set a\ = a[ and (72 = (T'2 in eq. (A. 8) and sum, 

^ {(Jia2s'3\siS2S'){aia2S3\siS2S) = Ss'^sJs'S • (5-12) 

a'i,(T2 

The sum over the third components of orbital angular momentum and S3 in the partial wave decomposition of 
eq. (5.5) becomes 

/ 2 7-1-1 

V {m's3f,\£'SJ){ms3f^\£SJ)Yir (p)rr(p)* = dt^e—— . (5.13) 
^ — ' 47r 

m',m,S3 

Here we have made use of the symmetry properties of the Clebsch-Gordan coefficients and of the addition theorem 
for the spherical harmonics [84], 

{m's3„\£'SJ) = (-1)^3+5 l^-^^j {-ss,^m'\SJ£') , 

y i^r(p)i^ = — • (5.14) 

m 

Whence, the sum of partial waves that matters for eq. (5.1) can be expressed as 

J2 T^lZ (P, a; ^) = (2-^ + ^' S)x{S£lf{a^a2h\hhlf , (5.15) 

CTi,CT2 i,j,e,s 
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with xiSil) defined in cq. (A. 7). Inserting the previous equation in eq. (5.9) the foUowing expression for Ilj results 

dTj,\i,e,s) 

(5.16) 



nr = -i^..,3 0e^'="''e^'=^''5:(2J+l)x(5a)2|Go(fcr),Go(fe); 



- (J0(«;i)n(j^0(«2)n 



Ily is an S-wave isovector self-energy contribution. This should be expected and it is due to the presence of the 
WT vertex for the coupling of the in- and out-going pions with a nuclcon, see diagram 7 of fig. 1 and cq. (5.4). 

We now consider the diagrams 8 in fig. 1, that involve the Born terms of pion-nucleon scattering. They are 
similar to the diagrams 6, though the nucleon self-energy is now due to the in-medium nucleon-nucleon interactions. 
Making use of eq. (3.30) and then integrating by parts, we have 



n« 



■ 2 2 

2/2 90^ 



(2^ 



(5.17) 



where the first term on the r.h.s. of the previous expression is isovector and the last one is isoscalar. The former 
is referred to as and the latter as H". Taking into account eq. (5.15) one is left with 



(2^)4 

dTJ,\£,£,S)\^ g\ q2 



^(2j+i)^(5a)2|Go(/ei)pC?o(fc2) ^ 



p — "^^J? — ^ ~ Go(A;i)„Go(fc2)n 



J,l,S 
j4 7,^ r j4 



+ x{SnfGa{ki)nGo{k2)nTj^{iJ,S) + Go{k^)pGo{k2)n x{StQfT%{e, I, S) + x{SnfT^j^{l, I, S) . (5.18) 



Eqs. (5.16) and (5.18) involve the knowledge of the derivative of the nucleon-nucleon partial wave amplitude with 
respect to the energy k^. Instead of the variable k\ we use the variable A, eq. (5.7), which is also the argument of 
Z/io and use the relation 



5/(a0) _ df{a^) _ ^df{a°) 



dk° 



dkl 



dA 



(5.19) 



with f{af) an arbitrary function that depends on and only through their sum. Let us now obtain an expression 
for the derivative of dTji/dA. For that, rewrite eq. (4.5) as 



Tji = Nji - Nji ■ Lio • Tji . 
Taking the derivative on both sides of the previous equation and isolating dTjj/dA, 



dTji _ 1 dNji_ _ 1 dNn , 
dA ~ " ■ dA " ' dA 



.o-D-jl-Njj-Djj.Njj. 



(5.20) 



(5.21) 



with 



D%{t, I S) = I + Nf,{£, I S) ■ L% , 



(5.22) 



the same matrix whose inverse is multiplying Nji(£,£, S) in cq. (4.31). Eq. (5.21) can be simplified by taking into 
account that Djj and Nji commute so that 



dT 



ji 



dA 



dNji 2 ^-^10 



(5.23) 
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At LO and NLO the previous expression reduces to 

dTjj 



dA 
dTji 



= D 



dA 



LO 



NLO 
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-N 



(0) 

JI 



2dL 



10 



dA 
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(0)" 

JI 
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(1)" 
JI 



dL 



(1) 
JI 



dA 



dL 



10 



dA 



D 



(1)" 
JI 



(5.24) 



with 



D 



(0) 

JI 



I + Nfi-L 



^(1) 
^Ji 



N^/i^)-Lw 



D 



(0) 

JI 



(5.25) 



Further, the standard notation {B, C} = B ■ C + C ■ B has been used in eq. (5.24). 

Eqs. (5.16) and (5.18) represent the contributions from diagrams 7 and 8 of fig. 1 to the pion self-energy in the 
nucleon medium. Their contributions are denoted by Ily and Hg, respectively. The former is purely isovector while 
the latter contains both an isovector and an isoscalar part, proportional to Sij^ and Sij, in that order, lly and llg'' 
are given by the same expression except by the global factor, proportional to for the former and to —g\<^l(t 
for the latter. This is just a consequence of the chiral expansion eq. (3.8) in the Born terms. On the other hand, 
is a N^LO contribution because it originates from the derivative with respect to k\ of the nucleon-propagator 
between the two pion lines. This propagator is not enhanced so that one order higher results as compared with 
the isovector part. 



6 Other nucleon-nucleon contributions and the cancellation of the 
isovector terms 

We now consider the calculation of those contributions that originate from the diagrams 9 and 10 of fig. 1, where a 

pion scatters inside a two- nucleon reducible loop. They are denoted by llg and Ilio, in order. As usual the diagram 
on the top corresponds to the direct part of the nucleon-nucleon scattering while that on the bottom represents the 
exchange part. The loop with the pions has to be corrected by initial (ISI) and final (FSI) state interactions, as 
denoted in the figure by the ellipsis which represent iterated nucleon-nucleon interactions. This iteration is the same 
as occurs for the nucleon-nucleon scattering in the nuclear medium, see fig. 6. The "elementary" nucleon-nucleon 
interaction Njj is dressed by the iterative process which gives rise to eq. (4.5), with Nji multiplied by the inverse 
of the matrix D,//. In this way, if we denote by £,ji{i,i, S) the elementary partial wave for a generic "production" 
process, Fjj{£,i, S), then the FSI dress it so that 

Fji{i,l S) = Dj}{i,l S) ■ ijiilX S) . (6.1) 

The matrix -Dj/, eq. (5.22), is already known from the study of the nucleon-nucleon interactions up to some order. 
On the other hand, ^j/ can be fixed following an analogous procedure to that used before for determining iVjj 
in section 4.1. In this way, ^j"-* is determined by expanding eq. (6.1) in powers of Lio up to (iio)" and then 
comparing with a full CHPT calculation up to ©(j/''^"). with at most n + 1 two-nuclcon reducible diagrams. Note 
that we have written m + n and n + 1 because for our present purposes the basic process, made up by a two- nucleon 
reducible loop with the two pions attached to one nucleon propagator, starts at 0{p~^), so that m = — 1, and it 
implies already one two-nuclcon reducible loop. In addition, both ISI and FSI are involved in the diagrams 9 and 
10 of fig. 1. Then, instead of eq. (6.1) we have 

Hji{i,i,s) = ^ D-l{i,i',s) ■ ^ji{e',e",s) ■ D-l{e",e,s) . (6.2) 

e',e" 

The LO result requires to employ Dj^j and to calculate the two-nucleon reducible loop to which the two pions 
are attached by factorizing on-shell the nucleon-nucleon scattering amplitudes. We use the notation D^'}-''*^ = 
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/ + Njj^'^'' ■ L5_o with n the chiral order, 



HM,o = dT"-^'j}-dT' ■ (6-3) 



Explicit expressions for DLiq are given below in eqs. (6.11) and (6.14). 



exact / / fact / exact / 7 exact / fact 



a) b) c) 

Figure 15: Diagrams that contribute to the calculation of ^^jj. Those two-nucleon reducible loops that contain the label 
"exact" must be calculated exactly in the EFT, while those with the label "fact" must be calculated with the on-shell 
factorization of the pertinent vertices. The filled circle in the figure indicates that the pion-nucleon scattering process 
contains both the WT and Born terms. 

At NLO one has an extra two-nucleon reducible loop. Expanding the Djj matrices in eq. (6.2) up to one iio 
and ^ji up to 0{p) we obtain 

eSn6V-2A^i?-iio-6r (6.4) 

Wc now match the previous equation with the result of fig. 15. In this figure wc have included inside each loop the 
labels "exact" or "fact" according to whether the loop is calculated exactly or by factorizing on-shell the nuclcon- 
nucleon vertices. The filled circle refers to the pion-nucleon scattering process that contains both the local and the 
Born terms, fig. 16. We denote by L^j) the two-nucleon reducible loop without external pions calculated exactly in 
CHPT and that occurs in figs. 15b and 15c. There is also the new contribution of fig. 15a whose exact calculation 
is denoted by DLj^j. The result is 

DL^j) - Ar(°; . DL,o ■ L« - l'JI ' DL,o ■ Nf^ . (6.5) 
The equality of eqs. (6.4) and (6.5), taking into account eq. (6.3) for implies that 

+ = DL^jI - {l'jI + • ^10, <^ } • DL,o . (6.6) 

In the last term we have the combination L^jj + (A^jj^)^ • iio which is 0{p) in our counting because it corresponds 
to the difference between an exact calculation of a two-nuclcon reducible loop and that obtained by factorizing the 
vertices on-shcll. The other contribution to S^^jj is given by DL'^jJ — ^jj*, as follows from cq. (6.6), that is also 0{p) 

by the same token. Finally, note that in the previous expression the two pions are attached to the loops DL^jj and 
DLio, while the remaining terms originate because of nucleon-nucleon scattering. 

\q /q q/\q 

k ^ J / A; ' ^ 



Figure 16: Born terms in pion-nucleon scattering. The vertices correspond to the lowest order pion-nucleon vertex. 
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The nucleon propagator before and after the fiUed circles in fig. 15 is the same so that it appears squared. This 
is required as the initial and final pion is also the same. We rewrite the nucleon propagator squared as 



0(C„-|pi-k|) 



+ 



^(IPi 



p° - fc? - E{pi - k)-ie^ p°-k°- E{pi - k) + ie 



d_ 
"d'z 



pl + z-k'i- E{pi -k) + i€ 



- + i{2-K)5{p\ + z-k°- E{p, - k))^(^„ - |pi - k|; 



(6.7) 



z=0 



The filled circles in fig. 15 consists of a WT pion- nucleon vertex and of the pion- nucleon scattering Born terms 
shown in fig. 16. Its sum is 



'2/ 



2/ ' 



E{pi - k + q) + ie 



kf - E{pi -k-q)+ie 



(6.8) 



We do not include the in- medium part of the nucleon propagator in the previous equation because for q° = 0{m^) 
the argument of the Dirac delta-function in eq. (3.2) is never satisfied as m,r ^ 0(nucleon kinetic energy). For 
the same reason, when performing the fcj-integration in the loop, the poles at fc? = Pi± — E{pi — k ± q), 
resulting from eq. (6.8), are not considered because the nucleon propagators will not be any longer of 0{p~^) but 
just of 0{p~^) (standard counting). A contribution two orders higher would then result. Once the A;"-intcgration 
is done the latter acquires from eq. (6.7) the value z+p\ — E{pi — k). The integration on ki for the evaluation of 
the two-nucleon reducible loop is analogous to the one performed in Appendix C for calculating the Lio function. 
The point is that Lio only depends on the energy of the external legs through the variable A = m{p1 + P2) — a^, 
cq. (5.7), that in turn only depends on the total energy. As a result, when the derivative with respect to z acts on 
a baryon propagator not entering in eq. (6.8), one has 



dz 



dpi 



m- 



dA 



dpi 



(6.9) 



Taking into account the chiral expansion of the nucleon propagators involved in eq. (6.8) for the pole terms and 
summing with the WT term, we have the leading contribution 



A 2 \ k 



(6.10) 



The antisymmetric tensor in eq. (6.10) gives rise to the isospin factor 2^3 in the evaluation of the loops in fig. 15. 
Notice that in the loop the pions are attached to the propagators of the two nucleons, the upper and the lower 
ones, and these two contributions sum symmetrically. We can then rewrite eq. (6.3) for this case as 



>(o) .mq° (. 2 , 
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W;iv 
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dA 



(6.11) 



In these equations we have included the subscript iv given their isovector character. In the same way for DL^j] 
one has 



DL 



(1) 

Jliiv 



.mq 
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2 9a-;;2 I ^3£ij3 



dL 



(i);i3 
JI 



dA 



(6.12) 



which corresponds to eq. (6.11) but substituting the term between brackets, where the nuclcon-nucleon vertices are 
on-shell, by its exact calculation. By applying eq. (6.6) we can fix C'jl^y in terms of eqs. (6.11) and (6.12). 

We now consider the case where the derivative with respect to z from eq. (6.7) acts on the baryon propagators 
involved in the Born terms of eq. (6.8) with ki =Pi + z — E{pi — k). The term E{pi — k) — E{pi — k ± q) can be 
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neglected when summed with for our calculation to NLO, so that the derivative with respect to z of eq. (6.8) 
yields the isoscalar contribution 

. (6.13) 

For any 13 the isospin identity operator gives rise to +2, instead of the factor 21^ of the isovector case. In this way, 
we can employ eqs. (6.11) and (6.12) substituting the vertex eq. (6.10) by eq. (6.13) and removing the action of 
the derivative md/dA. Thus, 



(0) _ 


gW 


JI;is 


Pil 




gW 


'10;is = 


P^l 


(1) _ 


gW 


JI;is 





Si, 



r>T — X r«3 

Here the subscript is is introduced given their isoscalar character. They give rise to IIjq. Notice that both DLiQ-is 
and DL^Jl^g arc one order higher than the analogous isovector terms in eqs. (6.11) and (6.12), respectively. This 
makes that starts to contribute at N^LO. 

Next, we proceed to obtain the expressions for the pion self-energy corresponding to the diagrams 9 and 10 of 
fig. 1 as a sum over partial waves. The leading contribution is obtained by using ^fj.j^y, eq. (6.11), with the result 

(n9+n-)|,,^ = -i!^^ {l-g\^)Ex{Siin2J+l) J J |^e-?V'=^''(Go(fcr),Go(fc.) 



•p 



where the part corresponding to IIj^q is the one proportional to g\ in the previous equation. A symmetry factor 
1/2 is included given the symmetry under the exchange of the two nucleonic external lines when they are finally 
closed. 

In ref. [1] we established that at 0{p^) all the contributions to the pion self-energy involving nucleon-nucleon 
interactions vanish {Vp = 2). This implies that the contributions IIt, n|", Eg and Il^g must cancel mutually at 
this order. Recall that the isoscalar ones, Ilg* and XIj^q, are 0{p^). The argument given in ref. [1] was a general 
one, without assuming any specific procedure for resumming the nucleon-nucleon interactions. We now show that 
UCHPT fulfills this requirement. When substituting into eqs. (5.16) and (5.18) the derivative dTji/dA eq. (5.24) 
at the lowest order, the following result is obtained 

(nr + n^^)L.o = ^^^(l-5if ) Ex(5a)^(2J + l) j^^j 0e-?V'=^''(^Go(fcr),Go(M. 

X • • ^ • [Dfr^r - G,{W,[k,UDf,'-r' ■ N'f ■ ^ • [^^?- V^) • (6.16) 

This equation is the same as eq. (6.15) but with opposite sign so that the cancellation with Ilg -|- II^q takes place. 

In the following of this section we work out several N^LO contributions that comprise the isoscalar term II^q as 
well those that are obtained by employing dT/dA to NLO and DL^j].^^ in Ily + Ilg' and Ilg + W-^q , respectively. For 
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the leading isoscalar contribution from IIio, which is ahcady 0{p ) due to the same reason as for ni*, we obtain 

■ • [DTr + xiseofiDTr' ■ Nff ■ [DTr'ij ■ (6-17) 

Including one more order in the calculation of Ilg, II^^^ and II^q requires the use of DL^jJ, eq. (6.6). The input 
functions DLj^j are given in eqs. (6.12) and (6.14) for the isovector and isoscalar cases, respectively. In this way, 



X 



(i);-ii-i 



'ar(l);-l o 



dA 



{L'r'+Nff.L-,lNf}}.'^\.[D^y'-^]-^\ . (6.18) 



8 W^LO' 



It is straightforward to show that fly+III^ calculated with dTji/dA evaluated at NLO with cq. (5.24), Ily + Ii; 
cancels with (llg + 11^^^) | ^^3^^, cq. (6.18). One has to replace Nj^j by its explicit expression in terms of L^jj , 

Nj^j = L^jj + Nj^j ■ Lio • Nj^j. Let us mention that there is an extra term for Ilg + 11™^ at 0{p^). It stems from 
a one more term to eq. (6.10) in the chiral expansion of the baryon propagators eq. (6.8). This contribution is 
suppressed by the inverse of the large nuclcon mass and will be considered when a full N^LO calculation of the 
pion self-energy in the miclcar medium is given. In Appendix B we evaluate explicitly DL^Jj and L^Jj needed for 
determining Nj^j and ^'j'j , eqs. (4.13) and (6.6), in order. Some steps introduced in the derivations of this and the 
previous section are also calculated explicitly. 

In summary, we have considered the calculation of the diagrams 9 and 10 of fig. 1, Ilg and IIio, in order. We 
have shown explicitly that up to and including NLO these contributions vanish with those evaluated in the previ- 
ous section, Ily and Ilg. Indeed, the cancellation obtained between these contributions in ref. [1] is more general 
because it was shown that once the Born terms in Eg and IIio are reduced to their leading contribution and summed 
to WT, as in eq. (6.10), the cancellation occurs. Some other contributions at N^LO have been also calculated, 
though they do not exhaust a full calculation to this order of the in-medium pion self-energy. It is then shown 
that at N^LO the just mentioned cancellation between the isovector contributions takes place as required by the 
findings of rcf. [1]. Note that at this order the full calculation of the two-nucleon reducible loops takes place. This 
clearly shows that the cancellation is beyond the factorization approximation valid at NLO, as should be according 
to ref. [1]. 



7 Nuclear matter energy density 

We study now the problem of the nuclear matter equation of state by applying eq. (2.5) up-to-and-including NLO. 
The diagrams required are shown in fig. 17. The first diagram corresponds to the energy of a free Fermi-sea. Its 
contribution, £1, is given by 

^l = Y^(Pp^p+Pna. (7.1) 

and is the only 0{p^) contribution. Nonetheless, since this is a recoil correction originating from the first term 
of eq. (3.5), one expects it to be suppressed numerically (it involves the inverse power of the hard scale m 
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Leading Order 




Vp = 1 

„6^ 



Next-to-Leading Order 




Vp = 2 



Next-to-Leading Order 




3.2 



Figure 17: Contributions to the nuclear matter energy up to NLO or 0{jP). The lines have the same meaning as in fig. 1. 



instead of ~ y/min^, the one in l/^ro-) In addition, as shown below, there is further a dimensional suppression. 
Hence, the NLO contributions involving the nucleon-nucleon interactions, suppressed by just one extra chiral order, 
could be of comparable size. This is of course an important remark for the possible saturation of nuclear matter. 

The NLO or Oij)^) contributions comprise the second and third diagrams. The former is the energy due to the 
one-pion loop nucleon self-energy whose expression, including a symmetry factor 1/2, is 



/ 



(27r)3 



4 /fco - £:(k) + it 



(7.2) 



The cut in fc° from the free one-pion loop nucleon self-energy is restricted to the lower half-plane of the A;°-complex 
plane, see eq. (3.16). In this way, by closing the first integral on the r.h.s. of the previous equation along an infinite 
semicircle centered at the origin and on the upper fc^-complex plane, the contribution from the free part of S^^ 
cancels. Concerning the last term, an analogous reasoning to that given previously in connection with fig. 3 can 
bo also applied here for the in- medium contribution S^^ ^ of S"^. This part is accounted for by the diagram 3.2 of 
fig. 17. Thus, the contribution that we keep now for the nuclear matter energy density from the second diagram of 
fig. 17, £2, is 



£2 = 1 



(27r)4 fco _ E{k)+ie 



^[0(Cp-|k|) + ^(e„-|k|)]S^. (7.3) 



Let us show that up to 0{p^) these two integrals give the same result. Taking into account the expression for the 
in-medium part of the one-pion loop nucleon self-energy, eq. (3.19), the first term on the r.h.s. of eq. (7.3) yields 
the integrals 



J W 
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-^(27r)5(A;° 



ml+iej (27r)4 fcO - E{]s.) + ie 
where the order of the integrations have been changed. We perform the shift k 



£;(k-q))0(e.3-|k-q|) , (7.4) 
k + q in the last integral, that is 
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finite. In this way, eq. (7.4) can be rewritten as 



/(0«..-W)(-.)/^ 



^q^-ml + ie -gO + £;(k) - E{k - q) + ze 



(7.5) 



If the higher order corrections from the difference -E(k) — E(\i — q) are neglected in the previous equation, the 
last integral is the one defining Ej, compare with eq. (3.16). Then, after summing over is, we have the same 
contribution as the last one in eq. (7.3) up-to-and-including NLO. So finally we can write 



d^k 
(27r)3 




(7.6) 



Let us evaluate now the diagrams 3.1 and 3.2 of fig. 17, that collectively are denoted as diagrams 3 in the 
following. These diagrams fully involve the nucleon-nucleon scattering. Their contribution to the energy density, 
£3, is 



(Tl ,(T2 ai ,Ct2 



In this expression we have explicitly shown the symmetry factor 1/2 and the sum over the spin (cTj) and isospin 
(ai) labels, with tj — >■ 0+. 

The LO in-medium nucleon-nucleon scattering amplitude calculated in section 4.2 does not depend on p^. The 

interaction kernel Njj only dc;p(?nds on at this order, while the resummation over the two-nuclcon intermediate 
states, that gives rise to Lio, Appendix C, depends on A and |a| as well. In the following we use as integration 
variables a and p introduced in eq. (5.6). The p°-integration from eq. (7.7) can be readily performed, with the 
result 



" dp" 

—Go{a + p)a^Go{a - p)a2 



^(|a + p|-^„je(|a-p|-^„J e{U-\^ + P\W{^o 



2a0 - E{ai + p)- E{a -p)+ie 2a° - E{aL + p)- E{3i -p)-ie 



(7.8) 



Here, wc have made use of the fact that only those terms with poles located in opposite halves of the p^-complcx 
plane survive after the p" integration. Those terms with the two poles located at the same half of the p^-complex 
plane vanish, as it is clear by closing the integration contour with a semicircle at infinite along the other half. We 
insert eq. (7.8) into eq. (7.7), and use the variable A, eq. (5.7), instead of aP. Thus, 



Cri,(T2 Oil,a2 



^3 2^ 2^ J (27r)3 (27r)3 27r ^"i"^^?.^'^) 

/Ti ITn n/1 rtli-, " \ / \ ' 



|a + p|-e«iM|a-p|-ea2) 
A - p2 -H ie 



g(^„, -|a + p|)g(^„, -|a-p|) 
A — p'^ — ie 



(7.9) 



where we made used that + /c^ — + a^)/2m. The resulting redefinition of rj is not indicated as it is not 
relevant for the following manipulations, as well as the factor exp(i?7a^) that is not shown. The first term between 
the square brackets on the r.h.s. of eq. (7.9) corresponds to the particle-particle part while the last one corresponds 
to the hole-hole part. Both of them originate by closing the nucleon lines in the diagrams 3 of fig. 17. Making use 
of 



^(|a + p| - ^„ J^(|a - p| - ^„ J =(1 - - |a + p|)) (l - ^(Ca 



.1)) 



(7.10) 
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Figure 18: Contours of integration Gi and Cjt on the yl-complex plane used to perform the integral in eq. (7.12). The 
former contour runs below the cut (dashed line) and the latter above it. The limits of the cut in A due to the hole-hole part 
of Lio, eq. (7.14), are indicated by A{a) and B{a). 



eq. (7.9) becomes 

^3 = -4^E E 



(27r)3(2^)3 27r'^ Ja,a.(P>a,A) 



(7\ ,(72 CKl ,(3:2 

^(eax-|a + p|) + ^(Ca, -|a-p|) 
Let us discuss the calculation of the integral 



A - p2 + ie 



- 2m5{A - v>'')e{U - |a + - |a - pL 



+00 



dA 



21: A - 



(7.11) 



(7.12) 



This integral is involved in the first two terms of eq. (7.11). As a preliminary result let us discuss the particle- 
particle and hole-hole parts in Lio, since they represent the inclusion of two-nucleon intermediate states in the 
medium. It is then convenient to express the function Liq employing the first form of the nucleon propagator in 
eq. (3.2). In this way 



■^10 



_(£k_ 
(2^ 

^(6 



Ia-k| 



6) 



£;(a - k) 



+ 



0(|a + k| 



-6) 



i;(a-k) + 



aP + - E{a. + \i)-ie aP + - E{ai + \i) + ie 



(7.13) 



Similarly as in eq. (7.8), only those contributions in eq. (7.13) with the two poles in lying on opposite halves 
of the fc'^-complex plane contribute. Then, 



-^10 



_fk_ 
(27r)3 



:|a - k| - ei)^(|a + k| - 6) ^(a - |a - k|)^(|6 - |a + k]; 



A--k^+i 



le 



y4-k2 - 



le 



(7.14) 



The first term is the particle-particle part and the last is the hole- hole one. Notice the different position of the 
cuts in A. While for the particle-particle case A has a negative imaginary part, —ie, for the hole-hole part the 
cut takes values with positive imaginary part, +ie. It is also worth mentioning that the extent of the cut in A 
for the holc-holc part is finite. This cut in the last term of cq. (7.14) requires k^ = A, but |k| is bounded so that 
the two 0-functions in the numerator are simultaneously satisfied. The extension of this cut is given in eq. (C.19). 
We denote its lower limit by A{a) and its upper one by B{a), corresponding to logarithmic branch points, where 
a = |a|. This observation is very useful for performing the integral of eq. (7.12). For its evaluation we consider 
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the contours C/ and C// of fig. 18. The dashed hne in the figure represents the cut in T^^^^J (p, a; A) due to that 
in the hole-hole part of Lio for A{a) < Re{A) < B{a) and positive imaginary part +ie. Physically it represents 
a real reshuffling of the occupied states by an in-medium pair of baryons respecting energy and three-momentum 
conservation. All the contours of integration include a semicircle at infinity centered at the origin along the upper 
part of the ^-complex plane. While the contour C/ runs along the real axis, and then below the cut, the contour 
Ci' runs above it, with the imaginary part +2ie. In both cases the pole in the denominator of eq. (7.12) at — ie is 
outside the contours of integration. Because of the convergent factor the integration along the semicircle at infinity 
is zero so that we can write 



+00 



dA 



JAr; 



277 ^ - p2 



:(p,a;A) 



Ci 



dA e^^^ 

p2 jg^aiS (P' ^) 



2tt A- 



Since the cut is outside the contour C/', see fig. 18, the integration along this contour is zero 

dA e^^-^ „ „ , 

T-„lll{v,a;A)=Q. 



2n A-p^ +ie 



(7.15) 



(7.16) 



Subtracting eq. (7.16) to eq. (7.15) we are left with 



dA 



c,2nA- p2 



-T^lZiP^a;A)- (f 

-L 



dA 



:(p,a;^) 



2tt A-p"^ + ie 

dA TglZ (p, a; A) - T^l^l (p, a; A + 2it) 



27r 



p^ +ie 



(7.17) 



Since the branch points are just of logarithmic type the integration along the vertical segments at A{a) and B{a) 
in fig. 18 do not contribute in the limit e ^ 0+. Notice as well that the limit 77 ^ 0+ is already taken in the 
last line of eq. (7.17). The amplitude T^l^l (P' a, A) can be obtained from the analytical extrapolation in A of the 
partial waves amplitudes 



(7.18) 



At LO Njj depends only on p^, although for higher orders it could depend also on 13, A and a^ in addition to p^. 
We can also apply here eq. (5.15), keeping explicitly the separation between the p^ and A variables, 



E 

f 1,<T2 



T^l^lip',a',A)^ J2 (2J+l)x(^«)'(ala2^3|/l/2/)'^}}(i?,^,5;p^a^A) . 
i,j,e,s 



(7.19) 



The analytical extrapolation in A docs not affect the expansion of the nucleon-nuclcon scattering in spherical 
harmonics associated to the angular variables which are left intact. From eq. (7.18) it follows that 



Ty,{p\ a^ A) - T}3 (p2, a^ A + 2ie) = N% + L%{b?, A) 



L\l{a\A) [ilUa^ ^ + 2te) LlMa^ A)] 



N% +L%{a\A + 2ie) 



VU^\A + 2ie) 



(7.20) 



In this equation we have taken into account that although Nji could depend on A for higher orders, in the difference 
A^}j(p2,a2, A) — Nj\{p^, a^, j4 + 2ie) this dependence cancels. The point is that the discontinuity in T}^ due to the 
right-hand cut is fully taken into account by multiplying the loop function Liq by the kernel Njj on-shcU. as in 
eqs. (4.5) and (4.31). The right-hand cut associated to the variable A is then removed in the process of calculating 
Nji order by order, as discussed in sections 4.1 and 4.2. 

As commented above, the difference L^^Q{a'^,A + 2ie) — V'^Q{a?,A) is due entirely to the hole-hole part of V^q, 
the last term on the r.h.s. of eq. (7.14). From eq. (7.14) we have 



L%{a\A + 2ie)-l}^^{a\A) = -m J ^^(^„, - |a + q|)^(e 

d\ 



q|) 



1 



i2TTm 



I 



{2ny 



^A-q^ + ie 
|a-q|)(5(A-q2) . 



A — q2 — ie 



(7.21) 
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Thanks to S{A — q^) the A-integration in eq. (7.17) is now trivial. "^^^ As a result eq. (7.11) turns into 
^3 = -4 E E C^J+mSil? I ^-|^^(C«.-|a + q|)e(^„,-|a-q|)(Ti3(q^a^q2) 



+ m 



d'p 1 - 0{U - |a + P|) - - |a - P| 



(27r) 



p2 _ q2 _ 



iV}3,(p^)->Ll3^(a^q^) . iV}3,(p^)->ii3(a^q2 + 2^e) 



(7.22) 



where we have indicated explicitly the integration variables in the different functions. The isospin index = ai+a2 
and for is = one should take just one the two possible cases with ai = —a2, \ai\ = 1/2. It is straightforward to 
show that ^3 given in cq. (7.22) is purely real, as it should be. First, the two ^-functions in the first line of eq. (7.22) 
imply that only the hole-hole part of Lio can have imaginary part. It follows that ilo(^^ + 2«e) = -^io(q^)* 
since, furthermore, Njj{p'^) + L^Q(q^) is a symmetric matrix (for the S = 1 and J = £±1 partial waves) or just 
a number (for the rest of partial waves), the diagonal elements of the product 



are positive real numbers. In this way we have for the imaginary part of £3, 



(7.23) 



Imr}3^(q2,a2,q2) 



+ mj A^{l-0(^„^-|a + p|)-0(C„,-|a-p|)}7^c5(p2-q2)T}3(p^a^q2)■^)Mp^a^q2)' 



(7.24) 



Taking into account eq. (4.31) and the expression for Lio, cq. (7.14), the imaginary part of TJ^ can also be calculated 
in terms of that of the hole-hole part of Liq. Substituting the result in eq. (7.24), it follows that 

Im(f3) = -4 E E (2-^ + mSU? J -f^-^eiU - |a + <l\)e{U - |a - q|) m 1 ^^^(p' " q') 



i,j,e,si3=-i 



{ 1 - 9iU - |a + P|) - e(€a. - |a - p|) + 9{U " |a + pMU " |a - p|)}t}^, • T 



(7.25) 



•,s) 



The quantity between curly brackets in the previous equation is 

[1 - e{U - |a + pD] [1 - - |a - pD] = ^?(|a + p| - - P| - U) . (7.26) 

But given |a| and |q| satisfying simultaneously the ^-functions in the first line, corresponding to two Fermi-sea 
insertions, it is not possible that they also satisfy simultaneously the two ^-functions in eq. (7.26), corresponding 
to the particle-particle part. Note that due to the Dirac (^-function in eq. (7.25) |q| = |p|. As as result, it is clear 
that Im(£3) = once eq. (7.26) is inserted in eq. (7.25). 
Writing explicitly 

Ty,{p\ a^, q2) . T}3 (p2, a^, q')* = E ^' P', a^, q2)T}3 (£', S; p^ a^, q^)*, (7.27) 



'^'^The values of that satisfy the two in-medium ^-functions in eq. (7.21) for a given |a| are comprised in the interval [j4(a), -B(a)], 
which is the domain of the A-integration in eq. (7.17). 
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Figure 19: Cyclic permutation of the free nucleon propagators with Fermi-sea insertions. The labels on the potential lines 
are included to appreciate the permutation of lines in the graph. As usual the thick lines in fig. 19 refer to the insertion of 
a Fermi-sea. 



we have for eq. (7.22) 

1 



i,j,e,si3=-i J \ J \ J 

d'p 1 - e{u - |a + p|) - eiu - |a - p 



^9{u-\^+<i\nu-\^-<i\) 



T}^,(f,f,5;q2,a2,q2) 



+ m 



p2 _ q2 _ 



. (7.28) 



The process followed from eq. (7.11) up to here can be schematically drawn as in fig. 19. The particle-particle 
part that results by closing the external lines in the diagrams 3 of fig. 17 is transferred to a reducible two-nuclcon 
loop entering the in-medium nuclcon-nucleon partial waves. This corresponds to the integration in p in eq. (7.28), 
which is sandwiched between two partial waves amplitudes, one of them complex conjugated. The integration on 
p^ in eq. (7.28) is linearly divergent because of the first integral on the second line. The product 



J2 T'fii^, S; p^ a^ q2)T}3^(^', i, S; p\ a^ q^)* 



(7.29) 



tends to a constant for p^ — > oo. Let us discuss how to regularize this integral, in the same way as already done 



for the calculation of the Lij functions. 



^ E|^r5|(£^5; p2, a2, q2)T]^,(f , ^, 5; p2, a^, q2)* 

(27r)3 
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^ T}^j{i,i', S; p2, a2, q2)T]3,(f , S; p^, a^, q^)* 



(7.30) 



Let us denote by N'j'j.^{si^,q^) ^ the limit for p^ — >• oo of iV}^j(p^, a^, q^) ^. At LO this limit is a constant for 
each partial wave. Then, we rewrite the previous integral as 

r,2 
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2^ 
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^0 



dp 



p2 _ q2 _ 



7V7^(p2,a^q^) +LS,(a^q2) Af}3^(p^ a^ q^) + (a^ q^ + 2ie) 
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^7/,cx>(a^q')" +Lro(a^q^ + 2^e) 



(7.31) 
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Now, since N'fj -)> N'fj ^ + 0(|p|"^) the integral 
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27r^ Jo — — ie 
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7V}^,(p^a^q2) >Ll|,(a^q^) iV}^,(p^ a^ q^) ' + L\lia' , + 2ie) 



(7.32) 



is convergent. The remaining integral in eq. (7.31) is expressed in terms of the function S'(q^), eq. (4.7), 
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(7.33) 



with 

T}^,_^(£^£,5;a^q2) = [iV]^,_^(a^q2)"' + LiUa^q^)l"' = lim r}^,(€', 5; p^ a^ q^) . (7.34) 
To simplify the notation let us define the symbols 

^Pi = E ^Ji (A 5; p^ a^ , q2)r}3^ (£' , 5; p2 , a^ q2 )* , 

^ooe = Y.Tfj^^{i,i',S-,B?,ct)Tf,^^{e,t,S-,B?,c^y . (7.35) 



The function ^(q^) depends on the same subtraction constant go already employed in the study of the vacuum 
nucleon-nucleon interactions. The final expression for £^ in eq. (7.28) is then 

^3=4 E E(2^+1M^^^)'/ 70(0^(Ca.-|a + q|)e(^„.-|a-q|)[-W^,5;q^a^q2) 
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(7.36) 



The sum of eqs. (7.1), (7.6) and (7.36) gives our result for the energy density, £, in nuclear matter at NLO, 

£ = £i+£2 + S3 ■ (7.37) 

We evaluate eq. (7.36) using the in-medium nucleon-nucleon partial waves determined at LO in section 4.2. The sum 
over partial waves shows good convergence already for maximum J = 4 and we sum up to J = 5. The results for 
the energy per nucleon, E/A = £/p, are shown in fig. 20 for symmetric nuclear matter, left panel, and for neutron 
matter, right panel. The inserted point with errors on the left panel of fig. 20 corresponds to the experimental 
values for the saturation of nuclear matter E/A = (—16 ± 1) MeV and p = (0.166 ± 0.018) fm"^ quoted in ref. [53]. 
The dotted line in the right panel is the result for neutron matter from the many-body calculation of the Urbana 
group [54]. It employs realistic nucleon-nucleon potentials and a fitted density dependent three- nucleon force in 
order to reproduce the experimental saturation point for nuclear matter. The rest of the curves, from top to bottom 
in both panels, correspond to the values of 50 = —0.25, —0.37 and —0.5 m^, in order. We observe that our curves for 
symmetric nuclear matter does have a minimum with a value in agreement with the experimental one, — 16± 1 MeV, 
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Figure 20: (Color online.) Energy per nucleon, £ / p, for symmetric nuclear matter, left panel, and for neutron matter, right 
panel. The (magenta) dotted line in the right panel is the result from the many-body calculation of ref. [54] using realistic 
nucleon-nucleon potentials. The rest of the lines from top to bottom correspond to different values of go = —0.25, —0.37 
and —0.5 m^, respectively. The point on the left panel is the empirical saturation one for nuclear matter [53]. 



for go — —0.30. However, the position is displaced towards too low values of ^ = ^„ = ~ 150 MeV, too small by 
a factor 1.7 compared with the value ^ ~ 266±10 MeV [53]. For the case of neutron matter, the curves are repulsive 
and are larger than the calculation of the Urbana group [54] above some density. It is clear from fig. 11 and 12 that 
we are employing for the calculation of eq. (7.36) nucleon-nucleon partial waves that do not reproduce closely the 
Nijmegen data in several partial waves, see figs. 11-13. As commented above, we are just employing the iteration 
of the one-pion exchange and two four-nucleon local vertices. One needs more elaborate nucleon-nucleon partial 
waves. Indeed, there are many mutual cancellations involved in the case of symmetric nuclear matter, between 
the purely kinetic energy term, fi, and between the S- and P- waves in ^3, with £2 negligible small. Nevertheless, 
we find rather encouraging that our curves in fig. 20 can reproduce the main trends oi £ / p both for symmetric 
nuclear and neutron matter despite they are obtained employing in-medium nucleon-nucleon amplitudes calculated 
only at LO. We already pointed out in section 4.1 that the one-pion exchange has a too large tensor force which is 
reduced by higher order counterterms (in the meson exchange approach this reduction is achieved by the exchange 
of p-mesons [83].) In ref. [21] this point is emphasized in its study of nuclear binding because a large tensor force 
leads to less binding energy. Indeed, the partial waves ^Si-^Di and ^Pq have large matrix elements of the one-pion 
exchange tensor operator [21] and these partial waves are not well reproduced in our study at LO. We show in 
fig. 21 the different contributions to £ / p for symmetric nuclear matter and go = —0.37 m^, that corresponds to 
the solid lines in fig. 20. Namely, £2/ p and £3/ p are given by the dashed, dotted and dot-dashed lines, with 

the sum corresponding to the full curve. As expected, the contribution from £2, eq. (7.6), is very small since the 
derivative of S/ with respect to fc° is 0{p^), as discussed in section 3.'^^'^ The other terms, £i and £3 have similar 
size, though the former is 0{p^) and the latter 0{p^). This is due to the fact that the kinetic energy term is a recoil 

[2\ 

correction stemming from i2^j(r, eq. (3.5), being suppressed numerically by the inverse of the large nucleon mass to. 
Notice that £3/ p scales like pC, which introduces, compared with £1/ p, the additional power of ^ times l/27r^ from 
the density and Air/p from mC ~ rn/go. Both contributions have the same order of magnitude as the resulting 
factor 2^/7rp ^ 1. Additionally, there is also an extra suppression of the kinetic term contribution because of the 

^l^NoticG that S/(fc") = for fc" = and E(k) is very small. 
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Figure 21: (Color online.) The partial contributions to £ / p, for symmetric nuclear matter and go = —0.37 m\ are 
shown. The (black) dashed, (magenta) dotted, (orange) dot-dashed and solid lines are for £\/p, £2/p, £'i/p and their sum, 
respectively. The point is the saturation one for nuclear matter [53]. 



dimensionality of space. The point is that Ei contains the integral of d|k||k|'* — $^/5, while the extra factor of 
density in £3 goes like so that a numerical factor 3/5 = 0.6 is suppressing 81. Then, the cancellation between 
the kinetic energy term and the one due to the nucleon-nucleon interactions is a consequence of keeping the natural 
size for the chiral counterterms, of similar size than the one-pion exchange as seen in section 4.1. Of course, the 
precise value resulting from such cancellations depends on as shown in fig. 20. Additionally, the presence of such 
cancellation enhanced this dependence. E.g. for the neutron matter case the kinetic energy dominates the energy 
per baryon and the dependence on go is smaller indeed. £3 depends implicitly on go through the nucleon-nucleon 
partial waves. Additionally, there is an explicit dependence from the first term on the last line of eq. (7.36), go^oo£- 
The implicit dependence is due to the truncated solution of eq. (4.16), as discussed in detail in section 4.1. The 
explicit one should be also related to this truncation given the close similarity between them. To make this clear let 
us notice that the partial wave — Tj^(^, ^, 5; q^, a^, q^) appearing in eq. (7.36) can also be written from eq. (4.31) 



as: 



T]^,(^,^,5;q^a^q^) = - Nj, \c^^) + L^{e? ,ci^) 



-1 



= -E^/(^'^'"^'q''^''q')^/(^''^""^'q''^''q')*H/ '(/',€, 5; q2)+Lro(a^q^)*<5^'£) ■ (7.38) 

From the previous equation the term proportional to iio is 

-Lro(a^q2)*^T}^,(^,f,5;q^a^q2)^}^,(f,^,5;q^a^q2)* = {-Rel^^ + ilml^^)^,, . (7.39) 

Then, a similar dependence on as that of go^oof results from eq. (7.39) as —goZiq^. The sum of both is 
—gQ{Yjqi — Soof)- Thus, as a higher order calculation should dismiss the dependence on —gg'Sq^, by analogy, we 
expect this to be the case also for go^oo^- 

Another way of considering our power counting in eq. (2.5) is to use it for correcting order by order nucleon- 
nucleon amplitudes determined in vacuum. In this way, one can use better nucleon-nucleon partial waves, e.g. 
calculated at higher orders in momentum, and use in-medium corrections (whose calculation is always more cum- 
bersome than diagrams in the vacuum) at lower orders. Another interesting issue left for further work is to explore 
the three- nucleon force influence on £. This requires to consider the calculation of the energy density in the nucleon 
medium one order higher or 0{p'^), since Vp — 3. 
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7.1 Some phenomenology 



In this section, we will give up the strict power counting scheme employed so far and try to analyse the possible 
effects of higher orders by some phenomenologically guided parameter fine-tuning. This will allow us to better 
understand the results obtained for nuclear and neutron matter in comparison to other recent studies, as e.g. in 
refs. [31,33,35]. 

As a first exercise, let us vary the parameter go in order to improve the description oi £ / p for the case of neutron 
matter, so that our results agree better with the dotted line in fig. 20 corresponding to the sophisticated many-body 
calculation of ref. [54]. The dashed line in fig. 22 is obtained employing = —0.62 mj. The so obtained fine-tuned 
curve is very close to the results of ref. [54], even up to rather high densities (the deviation is then less than 10%.) 
Note as well that this result is obtained with a value of go still of natural size, in the expected range around 
—0.55 m^. However, if we employ the same go for evaluating E for symmetric nuclear matter the resulting curve 
has the minimum at its right position, p ~ 0.16 fm~'^, but the value of the energy per baryon is around —42 MeV, 
which is an over-binding by a factor 2.5. 




Figure 22: (Color online.) £ / p for pure neutron matter. The (magenta) dotted line corresponds to the result of ref. [54]. 
The (black) dashed line is obtained from eq. (7.36) with go = —0.62 m^. The (red) solid line represents eq. (7.40) with 
go ~ —0.62 m% and go ~ —0.65 m%. See the text for further details. 



To further analyse the density dependence of nuclear matter, we rewrite the contribution from the NN interac- 
tions introducing by hand a new parameter, so that 

^3=4 E E + / 7^7^e(U-|a + q|)0(ea. -|a-q|)[-^}^,(£,^,5;q^a^q2) 



/ cPp f l-g(ea, -|a + p|)-g(ea, -|a-p|) ^ 1 ^1 
+ go^^,-mJ ^1 p2_q2_,, S,, - -Eo.,|J . (7.40) 

Here we have distinguished between go, that corresponds to the parameter go that appears explicitly in 83 because 
of the diverging integral in eq. (7.30), and go, on which £3 depends implicitly through the dependence on the 
nucleon-nucleon partial waves. The idea is to mock up higher order effects by varying independently go and go 
and exploit the phenomenological implications of such a procedure. While go affects nucleon-nucleon scattering in 
vacuum, as discussed at length in section 4.1, go affects only the nuclear matter equation of state. E.g. employing 
go = —0.67 m^, which implies a change of 7% with respect to go = —0.62 used above, and that we also keep 
here, the solid curve in fig. 22 is obtained. The agreement for p < 0.2 fm"^ between our results and ref. [54] is 
almost perfect. 
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Symmetric Nuclear Matter 
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Figure 23: (Color online.) £/p for symmetric nuclear matter. The two (red) solid lines correspond from top to bottom 
to (go, go) = (-0.977,-0.512) ml and (-0.967,-0.525) ml, in order. The (blue) dashed line is obtained from eq. (7.41) 
adjusting to the minimum position and value of the lowest of our curves. The (magenta) dotted line is the result of ref. [54]. 



For the case of symmetric nuclear matter the best results are obtained with go = —0.52 ml, and go = —0.97 mj. 
Again the magnitude of both numbers is of natural size. This is shown in fig. 23 by the two solid lines which 
have {go,9o) = (—0.977,-0.512) ml and (-0.967,-0.525) ml, in order from top to bottom in the figure. The 
corresponding value for E/A at the saturation point is —15.4 and —17.1 MeV, respectively. The experimental 
value given by the cross corresponds to — 16 ± 1 MeV. The position of the minima for the same lines is p = 0.169 
and p = 0.168 fm~^, compared to the empirical value p = 0.166 ± 0.019 fm~'^. In addition, the dotted line is the 
result from the many-body calculation of ref. [54] employing realistic nucleon-nucleon interactions, that includes a 
free parameter to fix the three-nucleon interaction in the nuclear medium to reproduce the saturation point. We 
observe that our results reproduce very well the saturation point and agree closely with the calculation of ref. [54] . 

As done in ref. [31] it is illustrative to compare our curves in fig. 23 with the following simple parameterization 
for the energy per baryon in symmetric nuclear matter 

- = (7-41) 
p 10m m'^ 

Interestingly, the nuclear matter incompressibility [85, 86] 

is correctly given once a and /? are known by adjusting the empirical nuclear matter saturation point. In this 
equation is the Fermi momentum at the saturation point. For the central values of the point in fig. 23, p = 
0.166 fm^^ and E/A — —16 MeV, the resulting nuclear matter incompressibility is K = 259 MeV, which is 
compatible with the experimental value K — 250 ±25 MeV [85]. Our curves can be also described rather accurately 
by eq. (7.41), as shown by the dashed curve in fig. 23 obtained by adjusting the minimum position and value of the 
lowest of solid curves. Both curves run very close to each other and start to deviate for densities above p ~ 0.25 fm~'^. 
The resulting nuclear matter incompressibility calculated from our results is K — 254 and 233 MeV, for the upper 
and lower solid curves, respectively. These values are compatible with the experimental value. The close agreement 
between our results and eq. (7.42) shows that, to a good approximation, the former admits an expansion in powers 
of the Fermi momentum as in cq. (7.41) for low f . 

In Table 1 we show the contributions to £3 (MeV) in nuclear matter for the different partial waves considered 
at the saturation point p = 0.16 fm^'^. The first column to the right of every partial wave corresponds to go — 
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-8.80 


-31.03 




-1.56 


-1.54 


^^1 


-42.85 


-38.39 


^G3 


0.47 


0.62 




9.37 


-1.42 




1.71 


1.76 


'D^ 


-9.31 


2.94 




0.60 


0.61 


'Pi 


10.20 


12.30 


'Fi 


-0.04 


-0.05 


'Pi 


2.48 


2.78 


'Gi 


0.29 


0.15 


'P2 


-1.13 


-1.13 




0.03 


0.028 


'D2 


0.07 


-1.65 




0.32 


0.35 




-0.28 


-0.28 




-0.38 


-0.38 




0.91 


1.11 


'h 


0.23 


0.26 



Table 1: Contributions to E/A (MeV) in nuclear matter for the different partial waves considered at p = 0.16 fm"^. 
The first cohinin to the right of every partial wave corresponds to go = —0.977 and the second one to go = 
—0.521 m^, with go fixed to —0.512 m^. The former is the top solid line in fig. 23 and the latter is nearly the 
dot-dashed line in the right panel of fig. 20. 



—0.977 and the second one to go = —0.521 m^, and go is fixed to —0.512 in the two cases. The former case is 
the top solid line in fig. 23 and the latter is nearly the dot-dashed line in the left panel of fig. 20. The kinetic energy 
contributes with 22.11 MeV per nuclcon. One observes clearly the dominant role of the S-, P- and ^Di waves. It 
is remarkable the large infiuence of the medium on the 'Po partial wave that despite being attractive in vacuum 
gives a repulsive contribution in the medium for {go, go) = (-0.977,-0.512) m^, and only slightly attractive for 
9o — go — —0.512 m^. The change of the ^Sq contribution with go is due to the fact that for go ~ —0.977 the 
phase shifts decrease with energy (they are 30 and 20 degrees at |p| ~ 100 and 300 MeV, respectively), instead of 
stabilizing at around 60 degrees like happens for go = —0.512 m^, see fig. 11. In this way, its contribution is less 
attractive. We also remark that the non-clastic partial waves, e.g. the 'Si^'Di, also contribute to ^3 from Ep£ 
and Soo^ and are included in the elastic ones. For a given partial wave they correspond to one term in the sum of 
two terms on £' in eq. (7.35). 

Thus, wc arc able to reproduce the equation of state of symmetric and neutron matter in terms of one fine-tuned 
free parameter, the value of go — —0.97 m^. for the symmetric nuclear matter case. The parameter go comes out 
always with a value close to the expected one, —mm.,^/ATT = —0.54 m^. For neutron matter one has go — —0.62 
and for symmetric nuclear matter go — —0.52 m^. In addition, it is worth stressing that the resulting values for go 
are all negative with the appropriate size, so it is justified to solve approximately eq. (4.16) making use of UCHPT. 
We think that this achievement is not a trivial fact and clearly indicates that our power counting, eq. (2.5), is 
able to map out properly the important dynamics in the nuclear medium and establish a useful hierarchy to allow 
for systematic calculations. At present, the strong dependence on ^iq of particularly the ^6*0 and 'Po partial wave 
contributions, see Table 1, seems to indicate that in order to go forward in the application of our chiral power 
counting eq. (2.5), one should consider the exact solution of Njj, eq. (4.16), and not its truncated solution in 
powers of g, eq. (4.10), as presently done. Notice that in comparison with refs. [31, 33, 35], that also applies Baryon 
CHPT to in-medium calculations, we employ the same approach both for vacuum and nuclear matter calculations 
and are able to compare with nucleon-nucleon scattering data. 

8 Conclusions 

In [1] wc derived a novel approach for an EFT in the nuclear medium based on a chiral power counting that 
combines both short-range and pion-mediated internucleon interactions. The power counting is bounded from 
below and at a given order it requires the calculation of a finite number of contributions, which typically implies 
the rcsummation of an infinite string of two-nucleon reducible diagrams with the leading multi-nuclcon CHPT 
amplitudes. These resummations arise because this power counting takes into account from the onset the presence 
of enhanced nucleon propagators and it can also be applied to multi-nucleon forces. However, ref. [1] did not 
discuss any specific non-pcrturbativc method for applying this novel counting scheme to practical calculations. 
We have developed in the present paper the required non-perturbative technique that allows us to perform these 
resummations both in scattering as well as in production processes. This non-perturbative method is based on 
Unitary CHPT, which is adapted now to the nuclear medium by implementing the power counting of rcf. [1]. 
We have first applied it to calculate the LO and NLO vacuum nucleon-nucleon interactions. There, we have 
also shown the tight connection between the UCHPT approach and dispersion relations for the nucleon-nucleon 
scattering amplitude. The former results as an approximate solution to the dispersive treatment of nucleon- 
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nucleon scattering in a chiral expansion of the imaginary part of the scattering amplitudes along the left-hand cut, 
taking advantage of the suppression of the two-nucleon unitarity loops along this cut. It was also shown that the 
subtraction constant go, employed for regularizing the reducible part of the two-nucleon reducible loops, realizes 
a convenient splitting between loops and trce-levcl contributions, while the exact solution for nucleon-nucleon 
scattering in vacuum has been shown not to depend on g^. For the in-medium case the LO nucleon-nucleon 
scattering is given as well. Then, the pion self-energy in nuclear matter up to 0{p^) was determined together with 
some other contributions at N^LO. The cancellation found in ref. [1] between all leading corrections to the linear 
density approximation for the pion self-energy is explicitly shown here for the amplitudes calculated utilizing the 
non-perturbative method developed here up to N^LO, which is a good check for the consistency of our approach. 
We have also addressed the calculation of the energy density of nuclear matter £. The non-perturbative technique 
developed gives rise to different contributions to £ whose imaginary parts cancel between each other and a real 
value results, as it should. We obtain saturation in symmetric nuclear matter and repulsion for neutron matter. 
The contributions from the nucleon-nucleon interactions are of similar size to those from the kinetic energy term, 
the latter being suppressed by the inverse of the large nucleon mass and a dimensional factor. It is remarkable 
that we obtain for go — — 0.62m^ a very good reproduction of sophisticated many-body calculations that employ 
realistic nucleon-nucleon potentials for the equation of state of neutron matter up to rather high nuclear densities. 
We can also achieve such a good agreement for the case of nuclear matter by allowing to distinguish between go 
and go, where the former parameter appears explicitly in the calculation of the energy per baryon while the latter 
appears implicitly through the nucleon-nucleon scattering amplitudes involved. By this splitting, we can mock up 
the; effects of higher orders. The parameter go can be fixed from the neutron matter equation of state. We then 
obtain a very accurate reproduction to £ as function of density in symmetric nuclear matter for go = —0.97 m^, 
with saturation at p = 0.17 fm~^ and E/A = —16 MeV, cf. the experimental values p = 0.166 ± 0.019 fm~^ 
and E/A = — 16 ± 1 MeV. Furthermore, the nuclear matter incompressibility comes out with a value between 
240-250 MeV, in perfect agreement with the experimental one of 250 ± 25 MeV. We interpret the success of our 
reproduction of the nuclear matter equation of state for both symmetric and neutron matter in terms of one free 
parameter for each of them as an indication that our power counting is a realistic one and that it is able to establish 
a useful hierarchy within the many contributions and complications inherent to nuclear dynamics. This opens the 
way to proceed systematically improving the calculations in a controlled way. 

Certainly, higher chiral orders should be worked out to obtain more precise results concerning the interesting 
problems considered here. In particular, it would be desirable to solve exactly eq. (4.15) in order to have free 
nucleon-nucleon partial wave amplitudes independent of go- In addition, a N^LO calculation will address the 
interesting question about the importance of three-nucleon forces for nuclear matter saturation. Similarly, a N^LO 
calculation of the pion self-energy is a very interesting task. It will merge important meson-baryon mechanisms like 
e.g. the Ericson-Ericson-Pauli rescattering effect [45] , with novel multi-nucleon contributions that can be worked 
out systematically within our EFT. More calculations and applications of the present theory to other interesting 
physical problems should be pursued. 
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A Partial wave decomposition of in-medium nucleon-nucleon ampli- 
tudes 



In this appendix, we derive the partial wave decomposition of the nucleon-nucleon scattering amplitudes in the 
CM frame. Our states are normalized as, 

1 — particle state: (p',i|p,i) = 5ij{2'K)^S{p' — p) 

2 - particle state: (p', jii2|p, iii2) = 5j,i,6j,i,{2n)^S{Pf - Pi) ^^"^ ^ 6{n - O') . (A.l) 

Here, Vf corresponds to the total four-momentum of the final state and Vi to that of the initial one, with W = 
V° = Vj, the total CM energy. In addition, Ei and E2 are the energies of the particles 1 and 2, in order. The 
indices i and j refer to any discrete quantum number required to characterize the states. In the CM frame the 
solid angle is denoted by fl and the modulus of the three-momentum by p = |p|. The two-particle states with 
well-defined orbital angular momentum are defined as, 

\£m,Hi2) = I (ipr/"(p)*|p,iii2) . (A.2) 

Taking into account eq. (A.l) it follows then 

{i'm' ,jij2\(.m,iii2) = ^^^^5e,'e,5rn'rri5J^^Jj^i^ . (A.3) 

The decomposition into states with well-defined total angular momentum J, third component fi, orbital angular 
momentum i and total spin S is given by, 

|p,c7ic72) = XI {'^i(^2S3\siS2S){mszn\lSJ)Yp{py\JiieSsiS2) , (A.4) 

where the indices cti and (T2 refer to the third components of spin, si and S2 to their maximum values and m to the 
third component of the orbital angular momentum l/^^^^ Next, we introduce the isospin indices ai, a2, for the third 
components, ri, T2, for their maximum values, and decompose the free state in terms of states with well-defined 
total isospin / and third component 13. In addition, the antisymmetric nature of a two-fermion (two-nucleon) state 
is introduced. 



(|p, cTiaicr2a2) - I - p, CT2a2fTiai)) = — ^ ^^^{(Jia2Si\siS2S){msiij\lSJ){a-ia2h\TiT2l) 

X y^™(p)*|J/X^5s3/i3) - ((72aiS3|S2Sl-S)(mS3MK5J)(a2aii3|T2Ti/)yr(-P)1-^M^^S3/i3)} , (A.5) 

with the repeated indices to be summed. This convention is used throughout this section. To simplify the no- 
tation we denote the left-hand-side of the previous equation as |p, o'iaiCT2Q!2)A, with the subscript A indicat- 
ing that the state is antisymmetrized. Applying the relations [84] p) = (— l)^Yg'"(p), {(T2(^is-?,\s2SiS) = 
(— l)'^~*^~*^(c7icr2S3|siS2<S'), and analogously for isospin, eq. (A.5) for the nucleon-nucleon case (si — S2 = ti = 
T2 = 1/2) simplifies to 

|p,aiaia2a2)A = \/4^ {<Ji<J2S^\siS2S){ms^ii\£SJ)Y[^{-p)*x{Sei)\JlitSssIh) , (A.6) 

7,5",.^, m, 1,13 

with 



72 lO l-\-S + I= even 



Strictly speaking in order to match with the normalization in eq. (4.6) we should multiply the right-hand-side of eq. (A.4) by 
^Jm/E. This is a relativistic correction of 0(p^). 
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In this way, x(5i?/) ensures the well known rule that S + i + I must be odd for any partial wave. Using the 
decomposition eq. (A. 6) we have for the scattering amplitude, 

a{p' , o-[a[a2a2\T{a)\p, aiaia2a2) A = 'iTT^^{(7[a!2s';^\siS2S'){m' s'^iJ,'\£' S' J'){aia2Ss\siS2S){ms3iJ,\£SJ) 

X (a;a'2Z3|TiT2/)(aia2Z3|riT2/)r;"'(p')i^r(p)*x('5'^'^)x('5^^)rj'j/(^'^V^) • (A.8) 

Here, Tjiji{£'S';£S) is the partial wave with final total angular momentum J', initial one J, final total spin 
S' , initial one S, isospin / and final and initial orbital angular momenta £' and £, respectively. Notice that in the 
previous equation we have distinguished between the final and initial total angular momenta J' and J, and similarly 
for the total spins S' and S. For free two nucleon scattering we have of course J' = J because of angular momentum 
conservation. This conservation law, the conservation of parity and the rule S + £ -\- 1 = odd imply that S' = S. 
However, the resulting matrix elements in the nuclear medium depend additionally on the total three-momentum 
of the two nucleons because the medium rest-frame does not coincide in general with their center-of-mass. This is 
why we have included a as an argument in the scattering operator, with the former defined in (5.6). Employing 
the orthogonality properties of the Clebsch-Gordan coefficients and spherical harmonics, one can invert eq. (A.8) 
with the result, 

^-Kx{S'£'I)x{S£I)Tj>ji{£'S';£S) = Y, J dp' J dp a(p', ala;a^4|T(a)|p, aiaia2a2)A(a'ia^4|siS25') 
X (m's^/x'|f5'J0(CTia2S3|siS25)(ms3M|^5J)(a;a'2«3|TiT2/)(aia2i3|TiT2/)y/(p')*y«"(p) • 

(A.9) 

This expression can be further reduced by making use of properties under rotational invariance so that the initial 
relative three-momentum p can be taken parallel to the z-axis. In deriving this simplification we omit the isospin 
indices that do not play any role in the following considerations, and introduce the symbol 

^A"S(P''P'") = A(p',cT;a;a^a'2|T(a)|p,aiaia2a2)A • (A.IO) 

Now consider the rotation i?(p), such that i?(p)z = p, that consists first of a rotation around the y-axis with an 
angle 9 and then a rotation around the z-axis with an angle (p, where 9 and (f> are the polar and azimuthal angles 
of p, respectively. We could also have taken first an arbitrary rotation of angle 7 around the z-axis. Then, 

i?(p)t|p,c7ia2;a) = D^\i^^ {R^)D^li^J {R'<)\pz,s^S2-,a'') , 

Sl,S2 

i?(p)t|p',a;a^;a) = J^^ D^/J^ {R^)D^/J? {R^)\p\-s[-s',; ^'') , (A.ll) 

with p" = i?(p)~^p' and a" = i?(p)~^a. The dependence on the total three-momentum has been made explicit in 
the state vectors to emphasize that the total three-momentum also is rotated. Inserting eq. (A.ll) into eq. (A.9) 

we have, 

Anx{S'£'I)x{S£I)Tj,jji£'S';£S) = Y^ J dp' J dpT!/J^^{p'',pz,a'')D<^/^f{R^rD'^J^V{^^^^ 

X y/(p')*yr(p)K'^2S3|siS25')(m'4M'l^'5'^')(<^i'72S3|siS25)(ms3MK5J) • (A.12) 
The spherical harmonics satisfy the following transformation properties under rotations, 

Y,7'ip') = J2D^UR^)Yfip"), 

m' 

YriP) = Y.D'im{B))Yt{i) . (A.13) 

m 

Inserting these equalities into eq. (A.12) wo arc then left with the following product of rotation matrices, 

D%V{R^rD%V{R^rDiUR^rDi'^'XR^)^^^^^^^^ (A.14) 
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We now take into account the Clebsch-Gordan composition of the rotation matrices [84], 

Y,D'-^L{R){m[m',M'\£,e2L)= ^ Dl^^l^{R)Dl^^l^{R){m,m2M\h£2L) . (A.15) 

M' mi,m2 

Since eq. (A. 14) appears in eq. (A. 12) times Clebsch-Gordan coefficients we can make use of the previous compo- 
sition repeatedly. First, 

^Di:/^)(i?t)D(i/2^)(i?t)(^,^2S3|siS25) =Y,Df,UR^){s,-S2a^\s,s^S) . (A.16) 

The rotation matrix D\, , , that appears on the right-hand-side of the first of the previous equalities, can then be 
combined in eq. (A. 12) such that 

(i?t)Z)^:)„,(i?t)(^'4^'|^'5'j') =Y^D^p^,{R^){fn'-a',n'\l'S'J') . (A.U) 

Similarly 

Y,Dill{R^)D^il{R^){mss,i\eSJ) = Y,Dl/J{R^){7nasll\eSJ) . (A. 18) 
Incorporating eqs. (A. 17) and (A. 18) in eq. (A. 12), the latter takes the form 

47rx(S'fI)x(Sn)Trjj{fS';£S) ^Y^J^P'/ ^^pT;/;/^(p^pz,a'0r/(p")*l^^(^)i?S(^^)4t^(^^) 
X (m'a^/i'|f 5V0(sls2^3|siS2.S0(™^3At|^5'J)(siS2^3|siS25) . (A.19) 

Let us first consider the vacuum case where the scattering amplitude does not depend on a. In this way the 
integration over p in the previous equation can be done explicitly taking into account the orthogonality relation 
between two rotation matrices [84]. For that let us recall our previous remark about the fact that an arbitrary 
initial rotation over the z-axis and angle 7 can also be included. In this way we take 

tr^^^ J dP^'j^'MR'TD^fiJiR') = j^S,,,S,,,Sjj, . (A.20) 
Inserting this back to eq. (A.19) one arrives at 

x{Se'I)x{S£I)Tji{e',e,S) = J dp''l^i|(p'',pz)r/(p'')(m'a^a3l£'5J)(0a3a3l^5J) 

X {s'iS2a3\siS2S){siS2a3\siS2S) . (A.21) 

In this expression we have made use that only m = gives a contribution to Yf^lz) and, as explained after eq. (A. 8), 
S' = S. In addition, we have also used that dp' = dp", since both vectors are related by a rotation. The subscript 
J' in Tjrjj is suppressed because J' = J and it is redundant. Also, we have employed the notation for the partial 
waves of section 4, Tji{i', i, S). 

We now come back to the in- medium case and keep the dependence on a. Here also to = so that p, = a^. Let 
us show first that a Fermi-sea with all the free three-momentum states filled up to ^ has total spin zero. This is 
required because for a given three-momentum pi one has two spin states that must be combined antisymmetrically 
because of the Fermi statistics so that 5 = for this pair. Then, since this happens for any pair, the total spin of 
the Fermi sea must be zero. Regarding total angular momentum we now give a non-relativistic argument to claim 
that the orbital angular momentum must also be zero. This is due to the fact that the nuclear medium in the CM 
of the two nucleons that scatter is seen with a velocity parallel to —a. In this way, both the CM position vector 
and the total three-momentum of the nuclear medium are also parallel so that their cross product vanishes. As a 
result, since the intrinsic orbital angular momentum of the medium is also zero, one expects that the total angular 
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momentum is zero for the system also in the CM frame of the two nucleons. Thus, J' = J also in this case and 
then, because of the same reasons as in vacuum, S' = S. Let us recall the remark after eq. (5.15) to justify that 
/ is conserved also in the nuclear medium. In addition the third component of total angular momentum must be 
conserved, fj, — ji' , and summing over /i one has 

A* 

given the unitary character of the rotation matrices. Then, 

x{Se'l)x{SiI)Tji{e',i,S) = ^^^^^^Yl J ^^"/ #"r^|(p",p5,a")(s'i44|siS25)(siS2a3|siS25) 
X Yp'{p")*{m'a'3a3\£'SJ){0a3a3\eSJ) . (A.23) 
This expression reduces to the one in the vacuum, eq. (A. 28), whenever the integral 

J2 J dp"T!r^,^{p",pz, a")(s'i44|siS25)(siS2^3|siS25)y/(p")*(m'4a3K'5J)(0a3a3KSJ) (A.24) 

does not dependent on a. In that case the integration over da" is simply 47r and eq. (A. 21) is restored. 

Eq. (A.23) can be further simplified because for the evaluation of a nucleon-nucleon partial wave amplitude 
one only needs to consider the direct term in the nucleon-nucleon scattering amplitude. This follows because the 
operator T is Bose-symmetric so that, 

^a7aS(p''P'a) = (p',ala;a^a'2|T(a)|p,(7iai<72a2) - (-p', a^aXa'i|r(a)|p, <7iaia2a2) . (A.25) 

When implementing the second or exchange term in eq. (A.23), re-including the isospin indices as well, and using 

the above referred symmetry properties of the Clebsch-Gordan coefficients and spherical harmonics, one is left with 
the same expression as for the direct term in eq. (A.25) except for the global sign —(—1)"^+^ '^^ . Summing both 
expressions the factor 

1 - (-1)^+^'+^ (A.26) 

arises. Given the definition of x{S£I) in eq. (A. 7) and imposing the rule that £ + S + 1 = odd and £' + S + 1 = odd, 
the factor x{S£I)x{Si' I) can be simplified on both sides of eq. (A.23). The latter then reads 

r/j{e',i,S) = ^(al<7^s'3|siS2S)(<7ia2S3|siS25)(0S3S3K5J)(m'4s3|^'5J) 

X {a[a2i3\TiT2T){aia2i3\TiT2T) j j dp' {p' ,a[a[a2a2\Td{a)\pz,aiaia2a2)YiV'-' {p')* , (A.27) 

with only the direct term, as indicated by the subscript d in the scattering operator. For the particular case of the 
vacuum nucleon-nucleon scattering the previous expression simplifies to 

Tjiie',i,S) = Y.^a[a'2s'3\siS2S){a^a2S3\s,S2S)i0ssSs\eSJ){m's'3S3\e'SJ) 
X {a[a'2i3\TiT2l){aia2i3\TiT2l) / dp' (p', a^aia^a'alTrflpz, c7iaiCT2a2)Y«r'(p')* • (A.28) 



B Explicit calculation of Eg and IIio up to 0{p^) 

In this section wo evaluate explicitly DL^^j and since they enter for fixing Nj^j and S^^Jj , eqs. (4.13) and (6.6), 
respectively. We also show with explicit calculations some of the steps introduced in the derivations of sections 5 
and 6. For the calculation of these diagrams one has to recall that the exchange of a wiggly line corresponds to local 
plus one-pion exchange terms, fig. 5. When the nucleon loop, to which the two pion lines are attached, includes 
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q q 

m m 






Figure 24: The two-nucleon reducible loop with only local vertices. The free part of the in-medium nucleon propagator in 
eq. (3.2) is indicated by a thin line while the in-medium part, proportional to the Dirac delta-function, is denoted by a thick 
line. 

only nucleon-nucleon local vertices, we then have T\q and Tn for the isovector and isoscalar cases, respectively. At 

the same order, when one of the nucleon-nucleon vertices in this loop corresponds to a onc-pion exchange then T\2 
and Ti3 result. Finally, when both vertices are due to one-pion exchange one has T14 and Tis. 

We first consider the two-nucleon reducible loop with only local vertices, fig. 24. For the isovector case the 
derivative with respect to z from eq. (6.7) acts on a nucleon propagator not entering in eq. (6.8). With the 
four-nucleon local vertex of eq. (4.2) we then have for Tig, 



Tio 



2/2 



{CS {Sa^ t^l'o) +Ct {^a^fa^afifiSm'oSi'; 



^ d f d^k ( 

-^a^,l3^ae,aSm'tSi'o) jT^og^ J (27r)'^ '^"^^^ ~ k)oGo{p2 + k)t X |Cs (5a„a^a</3^mofe 

where k = 1 — g\q^ /qq. The spin indices are indicated with Greek letters. The momentum integration is the same 
as for the function Lio, eq. (C.l), so that a factor dVi^/dA results. The previous equation contains both the direct 
and exchange terms, though the direct term is the one needed to evaluate the different partial waves. Notice that 
the previous equation only contributes to the S-waves. We now work out the spin projections of the direct term in 
eq. (B.l) with the result 



(B.2) 



As it should, we have then the combinations {Cs — 3Ct)^ and {Cs-\-Ct)'^ for 5 = and 1, respectively. The isospin 

projection corresponding to the operator ~l~ ^m'm 

rf,^) is 2^3, which excludes 7 = altogether. Keeping 

only the direct term in eq. (B.l) , we have 



*3 — 7^^^ 



3Ct) 



Fit 

dA 



(B.3) 



For the isoscalar contribution Tn ^ we just have to remove the derivative md/dA and replace the vertex of eq. (6.10) 
by eq. (6.13). The spin projection is the same as before, eq. (B.2). However, the isospin operator now is different 
and for the direct term just corresponds to twice the identity operator. As a result 



9W 
Pql 



{Cs + (45 



(B.4) 



and both S = and 1 contribute. In Tio and Tn there is no dependence on a, they depend only on |a|, so that 
one can calculate the projection into the S-waves as in vacuum, making use of eq. (A. 28). 

We now consider T12 and T13, with one local and one-pion exchange vertices. The isovector contribution 
corresponds to fig. 25, T/2, and fig. 26, TI2, depending on whether the one-pion exchange is between the final or 
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Figure 25: The internal four-momenta and discrete indices are indicated on the two diagrams. The figure on the left 
corresponds to t/j" and the one on the right to ■ 





Figure 26: The internal four- momenta and discrete indices are shown in the figure for TI2. 



initial two nucleons, respectively. For the sum of the diagrams in fig. 25 one has 



X {Tmm 

+ rff) Go{pi - k)mGo{j>2 + k)e^— — ^ , 



(B.5) 



whore the repeated indices are summed and r = p'l — pi + k. Instead of keeping Gq{pi — k)^^ we take 
— dGo{pi — k + ^)o/dz\^^Q, with A = (-^,0), as done in eq. (6.7). On the other hand, for the pion propagator 
we neglect its dependence on r^, since it is while is 0(p^). Then, the energy dependence enters in r/2 

similarly as in Liq and the derivative can be taken with respect to the variable A, cq. (6.9), with A ~> after 
the derivative is performed. Let us consider the isospin and spin projections for the direct term, given by the first 
square bracket in the previous equation. The isospin operator T^>rn'''t'i{'''mm + '^u) projection between states 
of well defined isospin given by 213. The spin operator, after some algebraic manipulation in the term proportional 
to Ct, reads 

) • (B.6) 

The structures ^a la^^a.iai and Ua ,am ' ^a.iat ai'e already projected in eq. (B.2) for the different spin states. 
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There is, however, the new structure (ctq^,q,„ • r)((?Q^/a< • r) whose matrix elements are 


S = , - , 


S=l 






/ 


-1 +1 \ 


ll-S^sall = 


-1 



I +1 


r| -V2rs{ri + ir2) (n + ir2Y 
-V^rs{ri - ir2) - 2r| \/2r3(ri + ir2) 

(ri - ir2)^ V^rsiri - 1x2) r| / 



(B.7) 



In this matrix the rows correspond to the final third component of the total spin, S3, and the columns to the initial 
one, S3. Then, we can write for the direct term with total spin S 

TliT = (I) ' iCs - SCr)^ J |^Go(pr - fc)™^?o(.. + % 

Tigris',, ss) = is^ 



+ ml 



/ \ 2 
, Kq Szjs I 9A 



2/ 



{Cs + Ct)- 



md f d^k 



dAj (27r)' 



Go(pi - k)mGo{p2 + k)e 



1 



9 I 9 S., S3 5 



(B.8) 



using eqs. (B.5) and (B.7). In the integrals of eqs. (B.8) it is convenient to perform the shift of the integration 
variable 

(B.9) 



which implies that 



Q 

Pi-k^ — - k , Q =pi+P2 
Q 

P2 + k^ — + k , 



r = p[ - pi + k ^ p' + k , p' 



P'l - P'2 



(B.IO) 



For the one-pion exchange between the initial nucleons, fig. 26, we proceed in the same way followed for above. 



- 12 

Performing the same transformation in the integration variable as in eq. (B.9), the expressions for t{2^~^'^ can be 
used for Tl'^~^'^ with the exchange of p' — )■ p. 

We now introduce the functions Lu, L^^ and J defined as 



^"« = V (^ 



1 



4 (k + r)2 + ml 

^k k' 



(27r)4 (k + r)2+m2 

rf^fc k^y 



GoiQ/2-k)mGoiQ/2 + k)e , 

Go(Q/2 - k)mGo{Q/2 + k)e = L^^a' + L{^r' 

Go{Q/2-k)mGoiQ/2 + k)e 



(27r)4 (k + r)2 + ml 



(B.ll) 



with a = Q/2. These integrals are further discussed in Appendix D. In terms of the functions introduced in 
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eq. (B.ll) one can write the different matrix elements of tI^^ as 

, md 



■ / \ 2 

^i2,d = ^3 ^2 J C^s - 3C7t) — (Lio - m^Ln) 



rr'/>'S=l/i i\ • "'"1 



Ta 
11 



+ 2a3P3 



md 



^■^"^^(1,-1) = -*3^^ (f|) {Cs + Ct)^ [ia, - ^a,fLl- + {p\ - ^p',f{Ll! + 2L\, + Ln) 



+ 2ip[ - ip'^){a^ - m2)(Lrr + in)' 



(B.12) 



T(ir'(S^, 1) = -^3^^ (f^) ' V2{Cs + Cr)^ [{a, + ^a,) {iff as + (L^r + i?iy3} 



+ 2(p; + ip'^){a, + ^a2)(L^^ + ^n) 
7^/2^r'(0,-l) = -7^/2w='(l,0), 



/,s=i 

12, d 



T^/2w='(-l,0) = -T/2r(0,l) 



12,. 



^(-i,-i) = r/,7^(i,i) 



(B.13) 



The isoscalar contribution originates by taking the derivative of the intermediate nucleon propagator in fig. 16. 
We have the same expression as for T/2, eq. (B.5), but removing the operator md/dA and with the replacement 



9l lqP 



" p 'nn' ^ j-2 ^2 • 

As a result, the isospin operator changes and for the direct term it is given now by 

mTe'eGo{Q/2 — k)mGQ{Q/2 + k)e . 



(B.14) 



(B.15) 



The projection for a state with 23 = ±1 is 2(4/ — 3)Go(Q/2 — k)±i/2Go{Q/2 + k)±i/2- For is = there is now a 
contribution from eq. (B.15) given by 



2(47 - 3)- [Go(Q/2 - fc)+i/2Go(Q/2 + fc)_i/2 + Go(Q/2 - fc)_i/2Go(Q/2 + k)+y2] ■ 



(B.16) 
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Hence, instead of eq. (B.8), wo arrive at 

.2 ^2 / „ . \ 2 



X \ [Go(Q/2 - fc)„Go(0/2 + k)i + GoiQ/2 - k)eGo{Q/2 + k)m] , 

X i [Go(Q/2 - fc)„Go(0/2 + fc), + Go(Q/2 - fc),Go(0/2 + fc)„] , (B.17) 

such that m = i = ±1/2 for = ±1 and m = +1/2, ^ = —1/2 for = 0. Similar expressions to eq. (B.13) can be 
written. E.g. for 5 = one has now 

= (4^ - 3)^^ (ly) ' {Cs - 3Gt)(Zio - m^Zn) • (B.18) 
Here, the tilde indicates the symmetric form 

Lf- = I {Lf-{m,i)+Lt^-ie,m)) , (B.19) 

with m and £ given in terms of 13 as explained above. The expressions for T^f^ are the same as those worked out 

for with the replacement p' p, as noted above for the T12 amphtudes. 

Let us consider the partial wave projection of T12 and T13. As discussed at the end of Appendix A, we can still 
use eq. (A. 28), valid in the vacuum, if the integral in eq. (A. 24) does not depend on a. For T/2 and r/g, with the 
one-pion exchange between the final nucleons, this is clearly the case because and only depend on a through 
its scalar product with p'. Thus, there is no angular dependence left on a once the integration over dp' is performed. 
For the case when the pion is exchanged between the initial nucleons, fig. 26, the resulting and Tl^ do not 
depend on the final three-momentum p'. In this way, the integration over dp' cannot remove the dependence on a, 
This also implies that this diagram only can contribute to partial waves with £' ~ 0, that is, "^Si and ^Di — )• 
However, as remarked above after eq. (B.IO), the exchange p' p transforms T/2, r/g into TI2, T^g and vice versa. 
In addition, one has to notice the symmetry between p and p' in eq. (A. 9) for the partial wave decomposition 
It is then clear that the same partial waves result for the diagrams of figs. 25 and 26 with the exchange I' I 
Thus, we can still use eq. (A. 28) but using the diagrams with the pion exchanged between the final nucleons. The 
elastic partial wave is exactly the same for both diagrams and is equal to — > evaluated as 

discussed. We denote the corresponding partial waves 

by TI^.jMA^S) and T^^.jj{i,lS). For T12, that requires 

7 = 1, only the partial wave 7i2;0i(0, 0, 0) = 2T[2 qi(0j Oi 0) is not zero. On the other hand, for T13 both isospin 

combinations occur. Then one has the partial waves 7i3;oi (0, 0, 0) = 27^3.q;^(0, 0, 0), 7i3;io(0,0, 1) = 27^3.;^q(0, 0, 1), 

'7i3;io(2,0, 1) = 7^3.10(2,0, 1) and 7i3;io(0,2, 1) = 7^3.io(2,0, 1). We have the following expressions for Tjj{£,i, S), 
omitting the subscripts 12 or 13, 

TjiieJA) = ^Ity {/ #'^/(p') [T|=Ho,o)(ooo|aj)(ooo|aj) 

+ (T|=i(+i,+i) + r|=i(-i,-i)) (oiiiflj)(oiilftj)] 

- y dp' (r,-i(p')i^r'(-i.o) + i"/(p')i^r'(+i,o)) (ooo|lij)(i - io|aj) 

- 1 dp' (y,-^(p')T|=i(o,+i) + y/(p')T|=i(o,-i)) (oii|aj)(ioi|nj) 

+ y {Y^\^)Tt\-^.^^)^yh^)'rt\-^^-^)){^m^m-m^j)\ ■ (B.20) 
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Let us recall that in order to apply the previous equations, the vector p must be taken along the z axis. Namely, 
for the partial wave projections, we take the reference frame with the axes 



P , 

p X a 

sin/3 ' 
P X (p X a) 
sin/3 



pctg^-acsec/3 . (B.21) 



PI \ ,•' P'l P2 \ ,/' P2 

m o PI - k m' ^ t p2 + (i 

" 3 t P 2 + fc "1 o P I - k p'l 

« m m' 

a) b) 



Figure 27: The internal four-momenta and discrete indices are indicated on the two figures whose sum (Ictcrmiiics Ti4,d. 
Note that the pion labels and four-momenta are exchanged for the initial and final states separately between the two figures. 



Let us move now to the evaluation of T14 and T15 where both vertices in the two-nucleon reducible loop at which 

the two pions are attached correspond to onc-pion exchange. As in the previous cases we start by calculating the 
isovector case. We restrict ourselves from the beginning to the direct contribution, corresponding to the diagrams 
in fig. 27, whose sum is 

4 „0, 



■14,d ■ 



9A 



If J 2/2 



V m'o ' 



_3 _a c I a i; o 



om'oo'i't'te 



e't'te'tt'm 



X i^aam ■ k)(CT^a; • k)- 



The diagonal matrix elements of the isospin operator 



-Go{pi - k)oGo{p2 + k)t . 



\ rn' o 



.c 3 a c 
om' 00 'i't 'U 



•T-^ 3 a c 



Go(Pi - k)oGo{p2 + k)t , 



present in eq. (B.22), are given between states with well defined isospin as 

2*3 Go(pi - k)^iGo{p2 + 1 . 
which is zero for 13 — 0. With respect to spin we can rewrite, 

(1? • r)a,,/3(l? • r)a^,a{ff ' k)^^^ (ct • k)^^^ = (r • k)'^6a,,ae5a^,a^ - [(r X k) • CTa^,aJ [(r X k) • CTc«„ 
+ i{r X k) • CTa^,a<(r • k)(5a^,a„ + i{r x k) • aa^,a^{r ■ k)Sa^,ai ■ 

The matrix elements of the spin operators Sa^,a^Sag,ae ^^nd • v) a< • v) between states with well defined 

total spin were already worked out in eqs. (B.2) and (B.7), respectively. We have now in addition the operator 



(B.22) 

(B.23) 
(B.24) 

(B.25) 



whose matrix elements are 



-i 


V +1 



-2v3 \/2(t;i + iv2) 
y/2{vi -iv2) V2{vi + W2) 

s/2{vx-iV'2) 2u3 / 



(B.26) 



(B.27) 
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and for S — 0. We follow in the previous matrix the same notation as in cq. (B.7). Taking into account this 
matrix and eqs. (B.2) and (B.7) one can determine the operator of eq. (B.25) between two-nucleon states with well 
defined third component of total spin S. The latter are inserted into eq. (B.22) together with the isospin factor 
2i3Go±i/2Go±i/2, eq. (B.24). As a result, the amplitudes Tf_f^ and Tf^]-{s'^,S3) are determined. 

The isoscalar contributions arise by taking the derivative of the intermediate nucleon propagator in fig. 16 with 
respect to z, as discussed in eq. (6.13). We have the same expression as for ri4,d, eq. (B.22), but removing the 
derivative md/dA and with the replacement of eq. (B.14). The resulting isospin operator is now given by 

2T'^,Xmrinr!,G^{Pi - k)oG^{p2 + k)t . (B.28) 
One can work out straightforwardly its diagonal matrix elements between states with definite isospin, 

2(9 - 87)^16-0(^1 - k)oGa{p2 + k)t + Go(pi - k)tGo{p2 + A;)o| , (B.29) 

with o + t = is. Then, instead of eq. (B.22) one has now, 

Ti5,d = i (ly) ^^(9 - 87) J {Go(pi - k)oGo{p2 + k)t + Go{pi - k)tGo{p2 + k)o} 

r)(o^aa„ • k)((7^a< • k)-^— 5-7^— 5- , (B.30) 

with o + t = as before. The tensor integrals required by eqs. (B.22) and (B.30) that involve one intermediate 
two-nucleon state with two one-pion exchanges are calculated in Appendix E. 

To determine the partial waves for T14 and T15 we have checked numerically that eq. (A. 24) does not depend on 
a at the level of one per mil. This is similar to the numerical accuracy to which the in-medium integrations have 
been calculated. As a result, we can use the same equation for partial wave projection as in vacuum, cq. (4.4), but 
now in terms of Ti^^d and Ti^^d- This reduces the calculational load because eq. (4.4) has not the integration over 
a as in eq. (A. 27). We denote by Ti4;ji{i' ,i, S) and Ti5,ji{i' ,i, S) the resulting partial waves for T14 and T15, in 
this order. 

Summing over the previous partial waves it follows that 

= Tii,d + Ti3,d + ri5,d . (B.31) 
Comparing with eq. (6.14) it is straightforward to determine L^J}- 



C The loop function L 
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The function Lio is given by 



no = t 



^-^ + 27r;i^(a - I? - k|)^(gV2 - fc^ - - k)) 



^ + 27ri^(a k|)^(gV2 + fcO - E(^ + k)) 



(C.l) 



This integration corresponds to the loop in fig. 7 with total in-coming four-momentum Q. In the following we 
define, 

1 . . Q 



a=-(pi+P2) 2 



(C.2) 



The different contributions to Lio are calculated according to the number of in-medium insertions in the nucleon 
propagators, eq. (3.2). The ^''-integration for the free part, iio,/, is performed by applying Cauchy's theorem, 



Lioj - 



(27r)3 QO-}^-^+ie 

^ ' ^ m m ' 



-m 



f d^k 
J (27r)3 



)3 k2 - A - k ' 



(C.3) 
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with A given in eq. (5.7). One has to keep in mind in the following the +ie prescription in the definition of A. 
In order to emphasize this, we will write explicitly the combination A + ie in many integrals, though the +ie is 
already contained in A according to eq. (5.7). The result in eq. (C.3) corresponds to eq. (4.9) that is regularized 
according to the dispersion relation eq. (4.7), 



^wj — 9o 



.m\fA 
An 



(C.4) 



For the one-medium insertion, Lio,m the ^^-integration is done by making use of the energy-conserving Dirac 
delta-function in the in-medium part of the nucleon propagator, eq. (3.2). We are then left with 

d^k ^(6-|k-a|) + 0(^2-|k + a|) 



Lio,m = m 



(27r)3 



k"^ - A-ie 



Let us concentrate on the evaluation of the integral, 

d^k e{^i - |k-a| 



^io,m(6,A|a|) = m 



f d^k g(gi 
J (27r)3 k2 



(C.5) 



(C.6) 



-— 7 "I £i — a/a arctanh ^ 



^ — VA arctanh-""^ ^ ^ 



log 



\^\+ii?-A 



I A ■ 4|a| -(|a|-a)2-A 

Here, we have taken into account that the Heaviside function in the numerator implies the conditions 

|a| > a , 



|k| e [|a| - ^1, |a| -l-^i] , cos 61 e 

|a| < a , 

|k| e [0,a - |a|] , cos 61 e [-1,1] 
|k| e [a - |a|,,^i |a|] , cos 61 e 



a- 



2|k||a| 



2|k||a| 



(C.7) 



Despite the separation between the cases |a| > ^ and |a| < a, both give rise to the same expression in eq. (C.6). 
In terms of the function ^io,m(Ci) |a|), eq. (C.5) reads 



iio,m(6,C2,A,|a|) =fio,m(6,A|a|) + Vm(6,A|a|) . 
For the case with two medium insertions 

—im\rA 



^10,1 



87r2 



■y"dk0(6-|kVA-a|)0(6-|kVA + a|) . 



(C.8) 



(C.9) 



Here wc assume that ^2 > Ci • If the opposite were true one can use the same expressions that we derive below but 
with the exchange ^ O ^2- This is clear after changing k — >• — k in the integral of eq. (C.9). Denoting by Q the 
angle between k and a the two step functions require 



2|a|\/A 



(C.IO) 



2|a|A/]4 

One has to impose that y\ < +1 and j/2 > —1, otherwise cos^ is out of the range [—1,-1-1]. In addition, it is also 
necessary that 2/2 > 2/i • 

yi<+l^|a|-6<^<|a|+Ci , 
2/2 > -1 ^ |a| - 6 < VI < |a| + 6 , 



yi<y2^ A< 



??+^2 ,„,2 



(C.ll) 
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The simultaneous consideration of the first and third of the previous conditions requires that (|a| — ^i)^ < A^ax 
for |a| > ^1, which in turn imphes 

|a| < . (C.12) 

Notice that the previous upper bound is larger than because ^2 > Ci- Prom eq. (C.12) it follows then that 
|a| — C2 < 0, and since it is always the case that (|a| + ^2)^ > ^max the second condition in eq. (C.ll) is satisfied. 
On the other hand, 

if |a| > ^ Ama. < (|a| + 6)' , 

if |a| < ^ Ar^a. > (|a| + 6)' • (C.13) 

For the final form of -Lio,d one also has to take into account the conditions, 

2/1 > -1 ^ \/I> ^1 - |a|, 

y2 < +1 ^ \/I > 6 - |a| . (C.14) 
Gathering together the conditions in eqs. (C.10)-(C.14) we have the following options, 

yi < -1 , 2/2 < +1 ^ 6 - |a| < VI < 6 - |a| , (C.15) 
which is not possible because ^2 > ^i- Also 

yi < -1 , 2/2 > +1 ^ vG4 < a - |a| . (C.16) 
This only holds for |a| < ^1. Then cos^ e [—1,+!] and -Lio,d = —im\/A/{2'K). Further 

-1 < 2/1 < +1 , 2/2 > +1 ^ - |a|| < \/I < min(a + |a|,C2 - |a|) . (C.17) 

In this case, cos^ G [2/i:+l] and Lio,d = —im{(( — {\/A. — |a|)^)/(87r|a|). It follows that ^1 + |a| < ^2 — |a| for 
|a| < (^2 - ^i)/2 and ^1 + |a| > ^2 - |a| for |a| > (^2 - ^i)/2. In both cases [min(^i + |a|,^2 - |a|)]^ < Amax, as can 
be easily seen. The last possibility is 

-1 < 2/1 < +1 , 2/2 < +1 ^ 6 - |a| < \/A < a + |a| . (C.18) 

For this case to hold, it is necessary that |a| > — But then Amax < (£,1 + |a|)^ so that the allowed upper 

limit for \/A is y/~A^^ not ^ + |a|. In this case, cos^ e [2/1,2/2] and Lio.d = —im{^i + ^1 ~ — 2|ap)/(87r|a|). 
In summary, 

v^<Ci-|a| , |a|<ei 
-8^(?? - (VA- |a|)2) , la - |a|| < VI < a + |a| , |a| < 

-8^(e?-(^/I-|a|n, \^i-M<^/A<^2-\b.\,^<\3.\<^ ' 
-8¥r(^? + ^2 - 2A- 2|a|2) , 6 - |a| < ^/A < y/j;;:^ , Sa^ < |a| < «i±«i 

As a technical detail one can take directly the derivative of eq. (C.19) with respect to A for a given interval. 
The different intervals do not imply the appearance of delta functions of A when differentiating as one check by 
evaluating directly the derivative of eq. (C.9). 



^10,d 



D Calculation of Lii, and L"? 

We now consider the integrals defined in eq. (B.ll), though we calculate here completely only the scalar integral 
ill. For the calculation of and L^^ one resorts to the Passarino-Veltmann reduction and explicit calculations 
and expressions are given in ref. [63]. Nevertheless, we evaluate at the end of this Appendix the free part of the 
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tensor integral L^i in order to explicitly show how to handle the rcgularization of a two-nucleon reducible loop 
different to g{A) in terms of 50- 

We evaluate Lu according to the number of in-medium insertions. After performing the A;°-integration by 
applying Cauchy's theorem we are left for the free part with 

, £k I 1 m f\ 1 

-Lll,/ = ~m I ,^ tt; ^r^:^ : — = — — / dx- 



(27r)3 (k + p)2 + m2 k2 - ^ - ie 8tt Jq [p-^x{l - x) + mix - A{1 - x) - le]^^^ 
^^■log -4-(|p|+»m.)2 ^^^^^ 



SttIpI ^rnl + {VA- |p|)2 ' 

In the previous equation, we have introduced a Feynman integration parameter x as an intermediate step and 
V— a ± if = ±z\/a for a > 0. For the one-medium insertion case, the fc'^-integration is straightforward due to the 
presence of the energy Dirac delta-function 

. ^ fj^ ^;n(a-k) / d^k e^ja-k) 

"'"^ 7 (2^)3(k2-A-ze)((k + p)2+m2) ^""y (27r)3(k2-^-ze)((k-p)2 + m2) • ^ 
Both terms in the sum can be obtained from the function 

/ M'fc g(6-|a-k|) 

«ii,m-my (2^)3(k2_^_i,)((k + p)2+^^2) • ^^■'^) 
Let us work out the scalar product k • p. For the integration, we introduce the reference frame 

z = a , 

a X p 



Prom the last relation we have. 



|a||p| sin /3 

y = z X X = actan^ — pcsec^ . (D.4) 



p = a cos ,6 — y sin /3 , 
k • p = |k| 1p|(cos cos /9 — sin^^ sin(/)sin /3) , (D.5) 



with and integration variables in eq. (D.3) and cos /3 = p • a . Let us perform the (/(-integration, 

1 



that also has been checked numerically. In eq. (D.6) we have 

a = k^ + p^ + 2|k||p|cos/3cos0 + m^ , 



(D.6) 



6= -2|k||p|sin/3sin6i , (D.7) 



where sin/3 = 1/I — cos^ fi. Next, we move to the cos^ integration of eq. (D.3). For this integration one has to take 

into account the presence of the Heaviside function, which implies the conditions already worked out in cq. (C.7). 
The latter determine an interval of integration for cos9 G [a;i(|k|), a;2(|k|)] for a given value of |k|. Then, we can 
write, 



nX2 J nX2 J 

/ dcos^— ^ = / dcosd—^ = 

J XI \/a^ — b^ Jxi va' + b' cos ( 



e + d cos2 e 

X2 



1 r, r6' + 2c'cOs6l „ r- ; TT^l 1 

^<^log ^ + 2Va' + &'cos6' + c'cos2(9 } , (D. 

vc' I L Vc' J J 
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with 



a' = (52-4kVsin2/3 , 
6' = 4|k||p|(5cos/3 , 
c' = 4k2p2 , 
(5 = + + . 

Now, we consider the final integration on |k| in eq. (D.2) and define the auxiUary function, 

X2(|k|) 



/ll,m(|k|) = 



Stt^IpI 

in terms of which eq. (D.3) reads 
|a| > a , 

4l,m(6,C0S^) 

|ai < a , 

4l,m(Cl,COS/3) 



log 



6' + 2c' cos e 



+ 1\Ja! + cosfii + c' cos2 6* 



xi(|k|) 



(D.9) 



(D.IO) 



|a|+6 



d|k| 



k2 - A - 



/ii,m(|k|) 



f Ki-lal /-lal+Cr 
\,/o ^6-|a| 



■d|k| 



k2 - - ie 



/ll,m(|k|) 



Then from eq. (D.2) one has 



ill,m = •^ll,m(^m,C0S/3) +4i,„j(^^,-COS/3) . 



(D.ll) 



(D.12) 



Here, we have only indicated those arguments that change for each term in the sum of eq. (D.2) to calculating 
iii.m- Indeed, l\\,m depends also on |a| and |p|. The integration over |k| when two medium insertions are represent 
can be done straightforwardly because of the additional Dirac delta- function of the energy that fixes |k| = \fA. 
Then, 



im^ r ,^(^(ei-|k\/A-a|)0(C2-|kVA + a|) 



t 

Sn^ J (k + p)^ + 



(D.13) 



We have the same ^-integration as in eq. (D.6), with the same result but now with |k| = \/A. The integration over 
cos 6* is the same as in cq. (D.8), though with a different integration interval for con 6 that is fixed by the values of 
A and |a|, according to the results of section C. They are collected here for ^ < ^2, 



cos^e [a;i,a;2] with [a;i,a;2] = < 



f [-1,1] , ^/A<a-|a| , |a|<a 

[2/1,1], |a-|a||<VI<a + |a| , |a|<i^ 

[yi, 1] , 16 - |a|| < \/l < 6 - |a| , < |a| < £i±^^ 

I [2/1,2/2] , 6 - ja] < \/3 < ^/A;^ , < ]a| < ^ 



(D.14) 



2 ^ i"i ^ 2 

with yi and 2/2 defined in eq. (C.ll). We also define, similarly as was done for £u^m, the auxiliary function 



/ll,d(6,6,COS/3) = g^2|p| 



b' + 2c' COS 6 I — 

^ V 2 V a' + V COS 6* + c' cos^ d 

yd 



(D.15) 



with |k| = ^/A for the values of a' , b' and d in eq. (D.9) and x\ and X2 according to cq. (D.14). In this way, 

iii,d(Ci,6,cos/3) = -«7r/ii,d(a,6,cos/3) . (D.16) 
For the case ^ > ^2 the change of variable k — )• — k is performed in eq. (D.13) and thus 

iii,d = -i7r/ii,d(6,6,-cos^) . (D.17) 
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Let us now consider the calculation of ^ corresponding to 

We now multiply the previous equality by S"^ and sum over repeated indices, then 

3^^!^, + L^lfP' = -"^/ (gs [(k + p)2 + -A-ie] ■ (^-^'^ 

Doing first the angular integration, one has 

11,/ 11,/^ 87r2|p|7o ik2-^-ze (|k| - |p|)2 + ^2 

Let us consider the expansion of the last factor in the previous equation for |k| — > oo, 

(|k| - |p|)2 + m2 |k| |k|3 
Adding and subtracting 4|p|/|k| to eq. (D.20) it results 

11,/ 11,/!- Stt^IpI Jo k^-A-te\ (|k| - |p|)2 + |k| 

The first integral is convergent because of the expansion in cq. (D.21) while the second one corresponds to g(j(A), 
eq. (4.7).#i2 

Next, we multiply eq. (D.18) by p"p^, sum over a and b and proceed in the same way as before. After performing 
the angular integration we have 

d^k (k • p)2 



)3 [(k + p)2 + m2][k2-^-ie] 
^ r^l1d '^'(^'+P'+"') ^^ k2 + p2 + m2 (|k| + |p|)2+rn4 \ 

As in eq. (D.21) we consider the limit of the last factor in the integrand for |k| — >■ oo 

_ k^ + p^ + ml (|k| + |p|)2+< _ _ V , om-^) (D 24) 

4|k||p| '°S(|k|-|p|)2 + m2 - 3k2+^^l^l ^"-^^^ 

Adding and subtracting the first term on the r.h.s. of the previous equation, similarly as done in eq. (D.22), we 
have a convergent and a divergent integral. The convergent one is given by 



m n k\k^ + p2 + ml) ( k^ + p2 + (|k| + |p|)2 + rn^ 4p2\ . mp\p^ + ml) 

Stt^Jo ' ' k^-A-ie V 4|k||p| °^(|k|-|p|)2 + m2 +3k2y * u^VA 



(D.25) 



while the divergent integral corresponds to 

-2 /■oo 1,2 _2 



The calculation of the integrals eqs. (D.19) and (D.23) explicitly shows how to regularize in terms of the 

subtraction constant go the linear divergences present in the calculation of the two-nucleon reducible loop integrals 
needed for the evaluation of the loops in figs. 25, 26 and 27. Since two pion exchanges are involved in the latter 
figure only the scalar integral with |k|'' in the numerator, denoted by ^^2^ diverges. Its divergent part can be 
straiglil f()r\\ai(ll> i<\u,ulariz('(l in Icrms of g{A), as in eqs. (D.22) and (D.26). Indeed, it just corresponds to g{A). 

#i2pj.Qj.jj practical point of view it is simpler to calculate algebraically the initial integral eq. (D.19) in cut-off regularization than 
the convergent one eq. (D.22). Then, the cut-off is sent to infinity while keeping only the divergent linear term. The cut-off can be 
expressed in terms of go by proceeding in the same way for g{A) in eq. (4.7). 
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E Calculation of L12 



The different integrals involved in the evaluation of T14 and T15 can be expressed in terms of a set of scalar integrals. 
The tensor structure of these integrals is determined by the matrix elements in eq. (B.25). We also perform here 
the shift of the integration variable k {pi —p2)/2 + k ~ p + k , as in cq. (B.9), and rewrite eqs. (B.22) and (B.30) 
accordingly. The integrals necessary for the calculation of the latter equations are evaluated in ref. [63]. Here, we 
only discuss the calculation of the basic scalar integral involved, L12. For integrals with more complicated tensor 
structure one uses the Passarino-Vcltmann reduction. In the expressions that follow it is always assumed that the 
A;°-integration has been done either by using Cauchy's theorem, for the free part, or employing the energy Dirac 
delta-functions from the in-medium insertions. On the other hand, since now two pion propagators are involved 
we join them in one introducing a Feynman parameter 



1 



dy- 



1 



with 



where 



[(k + p')2 + ml] [(k + p)2 + ml] Jo [(k + A)2 + M2] 



A = p' + (p - p')y , 
M^ = ml + 2y{l - y)(p2 - p . p') = + 2y(l - y){l - cos^)p2 



(E.l) 



(E.2) 



p • p = p cos if 



(E.3) 

In order to apply the results already derived in Appendix D, where only one pion propagator was involved, we take 
into account that 

^ -/t^ . (E.4) 

[(k + A)2 + M2]2 am2(k + A)2+M2 

as it follows from the definition of M"^ in eq. (E.2). The scalar function L12 is defined by 



^12 



d^k 



(27r)4 [(k + p')2+m2][(k + p)2 



em) 



]°/2-k° - E{ai-k)-ie Q° /2 - k° - £{& - k) 



e{^i-\a + k\) 



+ 



^(|a + k|-^^) 



_QV2 + fc°--E(a-Fk)-ie QO/2 + kO - Eia + ls.) +ie 
For the free part 



d 



^12,/ 



m 



dml 



dy 



"K ^0 



d^k 



(27r)3 [(k + A)2 + M2](k2 -A-ie) 



The integration over k was already done in eq. (D.l). Making use of this result one has 

1 



^12,/ 



m 



JO 



dy- 



[p^^ml-A- 2iM^7^ 



M 



(E.5) 



(E.6) 



(E.7) 



with M 



' M"^. For the part with one- medium insertion 

d^k ^™(a-k) 



^12, m 



+ m 



(27r)3 [(k + p')2 
r d^k 

J w 



■ ml] [(k + p)2 + ml] (k2 - A - ie) 
^7(a-k) 



)3 [(k - pO^ + m2][(k - p)2 + m2](k2 -A-ie) 



The two terms in the sum can be obtained from the function 



m 



dy 



d^k 



^™(a-k) 



(27r)3 (k2 - A - ie)[(k + A)2 + M2]2 



(E.8) 



(E.9) 
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This integral is analogous to £ii^m in eq. (D.3), but with p and replaced by A and M, in that order. Furthermore 
one of the factors in the denominator is squared. Following the calculation of in,m, we adopt the reference system 



z = a 



with 



a X A 

X = =; , 

|a||A| sin?7 
y = z X X = actanTj — A csecr? 



a • A = |p||a| [(1 — y) cos/3' + y cos/3] = |a||A| cosij , 
cos /3' = p' • a , 



(E.IO) 



(E.ll) 



so that the scalar product k • A = |k||A|(cos0cos?7 — sin 6* sin ^ sin ?^), where 6 and ^ are the polar and azimuthal 
angles of k. Let us perform the (^integration in eq. (E.9) 



1 Z"^'' 



1 



[(k + A)2 + M2]2 



1 /■^'^ 



1 



-1 



d 1 



2|A||k|cos6'acos77 27r 



2tt 



k^ + A^ + + 2|k||A|(cos0cosr7 — sin sin </> sin 77) 
1 



k^ + A^ + + 2|k||A|(cos0cosr? — sin sin sin 77) 



(E.12) 



The last integral is of the type already evaluated in eq. (D.6) where now 

a = 5 + 2|A||k| cos ^ cos 7/ , 
b = — 2|A||k| sin^sinTj , 

with (5 = k^ + A^ + = k^ + p2 + as in eq. (D.9). Then, eq. (E.12) reads 



2ir 



S + 2|A||k| cos 77 cos 



27r Jo [(k + A)2 + M2]2 



4A2k2 (cos2 9 - sin^ 77) + 4 1 A| | k| 5 cos 77 cos 6* + J 2 



3/2 



(E.13) 



(E.14) 



The cos 6* integration of the previous result includes the Heaviside function in eq. (E.9) that fixes the limits of 
integration to xi and X2 given by eq. (C.7). This integration is straightforward, see e.g. the similar one of 
eq. (D.8). Then, our result for ^12,™ is 



|a| > 6 
|a| < 6 

^12,m = 



d|k| 



|a|-«i 



k2 - A- 



7-/i2,™(|k|: 



te 



[Jo ^6-|a| J 



d\k\ 



k2-A- 



-/l2,m(|k|) . 



le 



(E.15) 



Here we have used the function 



TTllkl 



Tl2,m( 



dy 



2|A||k| COST? + 6 cos6 



(27r)2 Jo _ 4|A|2|k|2)^4A2k2(cos2 e - sin2 77) + 4|A||k|5cos77cos6i + S'^ 



X2(|k|) 



Xl(|k|) 



(E.16) 
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For the function Li2,m of eq. (E.8) wc have Li2,m = £i2.rn{£.m,cP, cf5') + ^i2,m(C£, -c/3, -c/3'), with c/3' = cos/3'. 

■im\/I /• - \■k^/A - a|)0(C2 - |kVA + a|) 



87r2 



[(k + p')2+m2][(k + p)2 



(E.17) 



The angular integrations are of the same type as aheady developed for the case of the one-medium insertion and, 
hence, we define 



/l2,d(VI) 



dy 



2|A|\/Acos?7-|-5cos^ 



(27r)2 Jq _ 4|A|2A)y^4A2A(cos2 6 - sin^ ,7) + 4|A|\/1^ cosr?cos6i + 52 



2(VA) 



i(VA) 



(E.18) 



with the integration limits given in eq. (D.14), where it is assumed that ^1 < ^2- In terms of this function 
L\2,d — —iT'fi2,d{^i,£,2, c/3, c/3'). When ^1 > ^2 we perform, as usual, the change of variables k — k in eq. (E.17) 
and then Li2,d = -i-^fi2,d{^2,^i, -0/3, -c/3'). 
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